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Abstract 

Macrogases are defined as two-component, large-scale celestial ob- 
jects where the subsystems interact only via gravitation. The macro- 
gas equation of state is formulated and compared to the van der Waals 
(VDW) equation of state for ordinary gases. By analogy, it is as- 
sumed that real macroisothermal curves in macrogases occur as real 
isothermal curves in ordinary gases, where a phase transition (vapour- 
liquid observed in ordinary gases and gas-stars assumed in macrogases) 
takes place along a horisontal line in the macrovolume-macropressure 
(O^fy^p) plane. The intersections between real and theoretical (de- 
duced from the equation of state) macroisothermal curves, makes two 
regions of equal surface as for ordinary gases obeying the VDW equa- 
tion of state. A numerical algorithm is developed for determining 
the following points of a selected theoretical macroisothermal curve 
on the (O^fy^p) plane: the three intersections with the related real 
macroisothermal curve, and the two extremum points (one maximum 
and one minimum). Different kinds of macrogases are studied in de- 
tail: UU, where U density profiles are flat, to be conceived as a simple 
guidance case; HH, where H density profiles obey the Hernquist (1990) 
law, which satisfactorily fits to observed spheroidal components of 
galaxies; HN/NH, where N density profiles obey the Navarro- Frenk- 
White (1995, 1996, 1997) law, which satisfactorily fits to simulated 
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nonbaryonic dark matter haloes. A different trend is shown by the- 
oretical macroisothermal curves on the (O-^v-^p) plane, according if 
density profiles are sufficiently mild (UU) or sufficiently steep (HH, 
HN/NH). In the former alternative, no critical macroisothermal curve 
exists, below or above which the trend is monotonic. In the latter 
alternative, a critical macroisothermal curve exists as shown by VDW 
gases, where the critical point may be defined as the horisontal inflex- 
ion point. In any case, by analogy with VDW gases, the first quadrant 
of the (O^Tv^p) plane may be divided into three parts, namely (i) the 
G region, where only gas exists; (ii) the S region, where only stars 
exist; (iii) the GS region, where both gas and stars exist. With re- 
gard to HH and HN/NH macrogases, an application is made to a 
subsample (N = 16) of elliptical galaxies extracted from larger sam- 
ples (N = 25, N = 48) of early type galaxies investigated within the 
SAURON project (Cappellari et al. 2006, 2007). Under the simpli- 
fying assumption of universal mass ratio of the two subsystems, m, 
different models characterized by different scaled truncation radii i.e. 
concentrations, Sj, Sj, are considered and the related position of sam- 
ple objects on the (O^Ty^p) plane is determined. Macrogases fitting 
to elliptical galaxies are expected to lie within the S region or slightly 
outside the boundary between the S and the GS region at most. Ac- 
cordingly, models where sample objects lie outside the S region and far 
from its boundary, or cannot be positioned on the (O^fv-^p) plane, are 
rejected. For each macrogas, twenty models are considered for differ- 
ent values of (Sj,Sj,m), namely Sj,Sj = 5, 10, 20, +oo (Sj,Sj, both 
either finite or infinite), and m = 10, 20. Acceptable models are (10, 
10, 20), (10, 20, 20), (20, 10, 20), (20, 20, 20), for HH macrogases, 
and (10, 5, 10), (10, 10, 20), (20, 10, 20), for HN/NH macrogases. 
Tipically, fast rotators are found to lie within the S region, while slow 
rotators are close (from both sides) to the boundary between the S 
and the GS region. The net effect of the uncertainty affecting ob- 
served quantities, on the position of sample objects on the (O^fv-^p) 
plane, is also investigated. Finally, a principle of corresponding states 
is formulated for macrogases with assigned density profiles and scaled 
truncation radii. 

keywords - galaxies: evolution - dark matter: haloes. 



1 Introduction 

Tidal interactions between neighbouring objects span across the whole ad- 
missible range of lenghts in nature: from, say, atoms to cluster of galaxies 
i.e. from micro to macrocosmos. The role of tidal interactions is of basic 
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importance in driving a wide variety of physical phenomena. In dealing with 
microscosmos, tidal forces between molecules are responsible for the occur- 
rence of the liquid and solid phase, and the presence of a triple point where 
the gas, liquid, and solid phase coexist, for an assigned homogeneous sub- 
stance (e.g., Landau & Lifchitz, 1967, Chaps. VII- VIII, hereafter quoted as 
LL67). In dealing with ordinary cosmos, the tidal action of a white dwarf 
star on a sufficiently close (filling the whole volume enclosed by the Roche 
equipotential surface) red giant companion, makes mass transfer into the 
white dwarf until a critical mass is attained and the star ends its life into a 
catastrophic Snla supernova explosion (e.g., Burrows, 2000). In dealing with 
macrocosmos, the tidal action induced by massive haloes on hosted galaxies 
affects their formation and evolution process, due to a larger depth of the 
potential well, resulting in a different correlation of observables with respect 
to galaxies in absence of massive halos (e.g., D'Onofrio et al., 2006). 

Ordinary fluids are collisional, which makes the stress tensor be isotropic 
and the velocity distribution obey the Maxwell law. Tidal interactions therein 
act between colliding particles (e.g., LL67, Chap. VII, §74). Astrophysical 
fluids (leaving aside extremely dense environments such as galactic nuclei) 
are collisionless, which makes the stress tensor be anisotropic and the veloc- 
ity distribution do not obey the Maxwell law. Tidal interactions therein act 
between a single particle and the system as a whole. 

Given that tidal interactions are at work in both collisional and collision- 
less fluids, the existence of an analogy between the two may be the subject 
of a legitimate question. To this respect, an investigation must necessarily 
be restricted to theoretical considerations, as astrophysical fluids (conceived 
as macrogases) cannot be tested in laboratory. More specifically, a macrogas 
equation of state has to be formulated in terms of three variables (macrovol- 
ume, macropressure, macrotemperature), and the related macroisothermal 
curves (i.e. the macropressure as a function of the macrovolume for selected 
constant macrotemperatures) has to be compared with their counterparts 
deduced from the van der Waals (hereafter quoted as VDW) equation of 
state for ordinary gases. If some analogy exists, it can be extended to (un- 
detectable) real macroisothermal curves as a working hypothesis. Finally, an 
application can be made to galaxies or clusters of galaxies. 

The VDW equation of state can be expressed in (dimensionless) reduced 
volume, reduced pressure, and reduced temperature. Similarly to the Lane- 
Emden equation for polytropes (e.g., Chandrasekhar 1939, Chap. IV, §4; 
Caimmi, 1986), the reduced VDW equation holds for a class of fluids instead 
of a single fluid. In general, the states of two systems with equal values of 
the reduced variables, are defined as corresponding states. According to the 
principle of corresponding states, two fluids which obey the reduced VDW 
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equation of state and exhibit equal values of two among three reduced vari- 
ables, necessarily exhibit equal values of the remaining reduced variable. For 
further details refer to classical textbooks (e.g., LL67, Chap. VIII, §85). 

In the light of an analogy between ordinary gases and macrogases, the 
formulation of a principle of corresponding states in the latter case could be 
highly rewarding. A macrogas equation of state was formulated in earlier 
attempts (Caimmi and Secco, 1990, hereafter quoted as CS90; Caimmi and 
Valentinuzzi, 2008, hereafter quoted as CV08), where the analogy between 
macrogases and VDW gases was only mentioned, and isofractional mass (m = 
const) curves were plotted for a few selected density profiles. The current 
paper aims to establish a closer analogy, where the macrovolume is related 
to the fractional radius, y, the macropressure to the fractional mass, m, and 
the macrotemperature to the fractional energy, <fi. The basic assumptions 
and the formalism remain unchanged with respect to the last parent paper 
(CV08). 

The present investigation is mainly devoted to the following points: (i) ex- 
pression of an equation of state for two-component astrophysical fluids, con- 
ceived as macrogases; (ii) comparison between macroisothermal curves and 
isothermal curves related to VDW gases, with regard to a simple guidance 
case and two cases which satisfactorily fit to observations or simulations; (iii) 
application to a subsample (N = 16) of elliptical galaxies (CV08), extracted 
from larger samples (N = 25, N = 48) of early-type galaxies investigated 
within the SAURON project (Cappellari et al., 2006, 2007, hereafter quoted 
as SIV, SX, respectively). 

The work is organized as follows. The equation of state of ideal and VDW 
gases are reviewed, and related isothermal curves are shown, in Section |2j A 
macrogas equation of state is formulated in terms of macrovolume, macro- 
pressure, macrotemperature, and related macroisothermal curves are shown 
for flat and steep density profiles, in Section [3j An application to elliptical 
galaxies for which masses, radii, and rms velocities can be determined, is per- 
formed in Section HI where the selection of acceptable models is made and an 
interpretation of the results is outlined. The conclusion is drawn in Section 
EJ Further details on two specific points are reported in the Appendix. 

2 Ordinary fluids 

Let ordinary fluids be conceived as fluids where the effects of gravitation on 
the equation of state may safely be neglected e.g., on the surface of the Earth. 
The simplest description is provided by the theory of ideal gas. 

Ideal gases are collisional fluids defined by the following properties: (i) 
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particles are identical spheres; (ii) the number of particles is extremely large; 
(iii) the motion of particles is random; (iv) collisions between particles or with 
the wall of the box are perfectly elastic; (v) interactions between particles or 
with the wall of the box are null. 

The equation of state of ideal gases may be written under the form (e.g., 
LL67, Chap. IV, §42): 

P V = kNT ; (1) 

where p is the pressure, V the volume, T the temperature, N the particle 
number, and k the Boltzmann constant. The product, pV, has the dimen- 
sions of an energy and, in fact, the mean kinetic energy per degree of freedom 
equals the product, (l/2)kNT, and the mean kinetic energy of motions along 
the X p axis reads: 

1„ o 1 



(E kin ) pp = -Nma 2 pp = -kNT ; (2a) 
fna pp = kT ; (2b) 

where m is the mean particle mass and a pp the rms velocity component along 
the X p axis. In the light of the theory of ideal gases, Eqs. (JTJ) and (j2J) disclose 
the meaning of the Boltzmann constant: for fixed pressure, volume, and 
particle number, the mean kinetic energy remains unchanged, regardless of 
the nature of the gas. 

In getting a better description of real gases, the above assumption (v) is 
relaxed and interactions between particles are taken into consideration. The 
VDW generalization of the equation of state of ideal gases, Eq. (pQ), reads 
(van der Waals, 1873): 

(p + A^\(V-NB) = kNT ; (3) 

where A and B are constants which depend on the nature of the particles. 
More specifically, the presence of an attractive interaction between particles 
reduces both the force and the frequency of particle-wall collisions: the net 
effect is a reduction of the pressure, proportional to the square numerical 
density, expressed as A(N/V) 2 . On the other hand, the whole volume of the 
box, V, is not accessible to particles, in that they are conceived as identical 
spheres: the free volume within the box is V — NB, where B is the volume 
of a single sphere. For further details refer to specific textbooks (e.g., LL67, 
Chap. VII, §74). 

The isothermal (T = const) curves for ideal gases are hyperbolas with 
axes, p = ^fV, conformly to Eq. (pQ). In VDW theory of real gases, the 
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Figure 1: Isothermal curves related to ideal (left panel) and VDW (right 
panel) gases, respectively. Isothermal curves (from bottom to top) corre- 
spond to T/T c =0.85, 0.90, 0.95, 1.00, 1.05, 1.10. No extremum point exists 
above the critical isothermal curve, T/T c = 1. 



isothermal curves exhibit two extremum points, which reduce to a single ho- 
risontal inflexion point when a critical temperature is attained, as shown in 
Fig.[TJ Well above the critical isothermal curve, T ^> T c , the trends exhib- 
ited by ideal and VDW gases look very similar. Below the critical isothermal 
curve, T < T c , the behaviour of VDW gases is different with respect to ideal 
gases and, in addition, the related isothermal curves provide a wrong de- 
scription within a specific region where saturated vapour and liquid phases 
coexist. Further details are shown in Fig.[2J Above the critical isothermal 
curve (T = T c ) the trend is similar with respect to ideal gases. Below the 
critical isothermal curve and on the right of the dashed curve, the supersatu- 
rated vapour still behaves as an ideal gas. Below the critical isothermal curve 
and on the left of the dashed curve, the liquid shows little change in volume 
as the pressure rises. Within the bell-shaped area bounded by the dashed 
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curve, the liquid phase is in equilibrium with the saturated vapour phase. 
A reduced volume implies smaller saturated vapour fraction and larger liq- 
uid fraction at constant pressure, and vice versa. The VDW equation of 
state is no longer valid in this region. The dashed curve (including the cen- 
tral branch) is the locus of intersections between VDW and real isothermal 
curves, the latter being related to constant pressure where liquid and vapour 
phases coexist. The dotted curve is the locus of VDW isothermal extremum 
points. 

A specific (T/T c = 0.85) VDW and corresponding real isothermal curve, 
are represented in Fig. [3j The VDW isothermal curve and the real isothermal 
curve coincide within the range, V < Va and V > Ve. The VDW isother- 
mal curve exhibits two extremum points: a minimum, B, and a maximum, 
D, while the real isothermal curve is flat, within the range, Va < V < Ve- 
Configurations related to the VDW isothermal curve within the range, Va < 
V < Vb (due to tension forces acting on the particles yielding superheated 
liquid), and Vd < V < Ve (due to the occurrence of undercooled vapour), 
may be obtained under special conditions, while configurations within the 
range, Vb < V < Vd, are always unstable. The volumes, Va and Ve, cor- 
respond to the maximum value in presence of the sole liquid phase and the 
minimum value in presence of the sole vapour phase, respectively. 

The surfaces, ABC and CDE, are equal, as first inferred by Maxwell (e.g., 
Rostagni, 1957, Chap. XII, §19). The VDW and real isothermal curves 
represented in Fig. [3] being related to the same temperature, T, the cycle, 
ABCDECA, is completely both isothermal and reversible, and the work, W, 
performed therein cannot be positive to avoid violation of the second law of 
the thermodynamics. The cycles, ABCA and CDEC, occurring in counter- 
clockwise and clockwise sense, respectively, are also completely both isother- 
mal and reversible. Accordingly, IVabcdeca = Wabca — Wcdec < 0. A sim- 
ilar procedure, related to the reversed cycle, ACEDCBA, yields Wacedcba = 
Wcedc - IVcbac < 0. Then IVabcdeca = Wacedcba = °> which implies 
W^abca = W-cdec = H^cedc = W^cbac and, in turn, the equality between the 
related surfaces. For further details refer to specific textbooks (e.g., LL67, 
Chap. VIII, §85). 

In order to simplify both notation and calculations, it is convenient to 
deal with (dimensionless) reduced variables (e.g., Rostagni, 1957, Chap. XII, 
§16; LL67, Chap. VIII, §85). To this aim, the first step is the knowledge of the 
parameters related to the critical point, V c ,p c ,T c . Using the VDW equation 
of state, Eq. flSJ), the pressure and its first and second partial derivatives, with 
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respect to the volume, read: 




kNT N 2 

P= V^NB- A V^ ; iV = const ! W 
kNT A N2 

~~~{V-NBf +2A V^ 5 (5) 

2kNT N 2 
yep = {V - NBf - 6 V 5 (6) 

where the domain is V > NB, V = NB is a vertical asymptote, and p = 
is a horisontal asymptote. The critical isothermal corresponds to the highest 
temperature allowing a liquid phase, which occurs therein only at the crit- 
ical point. The critical isothermal curve exhibits neither a minimum nor a 
maximum, which are replaced by a horisontal inflexion point coinciding with 
the critical point. Accordingly, (9p/dV)v c ,T c = 0, (d 2 p/dV 2 )v c ,T c = 0, and 
p c = kNT c /(V c - NB) - AN 2 /V 2 . The solution of the related system is: 

V C = 3NB ; (7) 
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Pc = 7^—2 : ( 9 ) 

(10) 



A;A^T C 8 ' 

where, in general, the compressibility factor, Z = pV/(kNT), defines the 
degree of departure from the behaviour of ideal gases, for which Z = 1, 
according to Eq. (JTJ. For further details refer to specific textbooks (e.g., 
Rostagni, 1957, Chap. XII, §20; LL67, Chap. VIII, §85). 
With regard to the reduced variables: 

V v T 

the ideal gas equation of state, Eq. (pQ), and the VDW equation of state, 
Eq. reduce to: 

(12) 

= lr ; Y>\; (is) 
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and Eqs. (SJ) , (JSD , and (jSJ), reduce to: 

8 ^ 3 
^~ 3f - 1 ~ ~f 
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(14) 

24 f 6 

+ ^ ; (15) 



V^r; r , r (3f-i) 2 p 

144 y 18 
(3f"l) 3 f 4 ' 1 J 

where, for assigned JT, the domain of the function, tfCf), is f > 1/3, = 
1/3 is a vertical asymptote, and $ = is a horisontal asymptote. In the 
special case of the critical point, f = 1, = 1, $f = 1, the partial derivatives 
are null, as expected. 

The extremum points, via Eq. (1151) . are defined by the relation: 



f^f) = %^ l = r; (i7) 

which is satisfied at the critical point, as expected. The function on the 
left-hand side of Eq. ( ITT)) has two extremum points: a minimum at f = 1/3 
(outside the physical domain) and a maximum at f = 1, where JT = 1. 
Accordingly, Eq. flTTJ is never satisfied for f > 1, which implies no extremum 
point for related isothermal curves, as expected. The contrary holds for 
JT < 1, where it can be seen that the third-degree equation associated to 
Eq. (TT5|) has three real solutions, related to extremum points. One lies outside 
the physical domain, which implies f < 1/3. The remaining two are obtained 
as the intersections between the curve, f(f), expressed by Eq. (ITT)) , and the 
straight line, y = ft, keeping in mind that /(1/3) = 0, /(l) = 1, and 
lim f _ +oo /(f) = 0. 

The third-degree equation associated to Eq. ( TT5)) . may be ordered as: 



f 3 -9af 2 + Qaf -a = ; (18a) 

a = ; (18b) 

with regard to the standard formulation (e.g., Spiegel, 1968, Chap. 9): 

x 3 + aix 2 + a 2 x + a 3 = ; (19) 
the discriminants of Eq. (1 18 aft are: 

„ Sao — On , , 

Q = Q 1 =a(2-9a) ; (20) 
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9aia 2 - 27a 3 - 2a? a(l - 18a + 54a 2 
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D = Q 3 + R 2 



[1 - 4a) 
~~ 4 



(21) 
(22) 



where D = in the special case of the critical isothermal curve ( JT = 1, a = 
1/4), D < for JT < 1, and D > for > 1. Accordingly, three coincident 
real solutions exist if D = 0, three (at least two) different real solutions if 
D < 0, one real (outside the physical domain) and two complex coniugate if 
D > 0. 

The three real solutions (D < 0) may be expressed as (e.g., Spiegel, 1968, 
Chap. 9): 



(23a) 
(23b) 
(23c) 




3 Gl 



-ai 



(23d) 



where a\ = 
a = 1/4, Q 
fc = % = 



—9a and, in the special case of the critical isothermal curve, 
= -1/16, D = 0, which implies % = min(f 1 , f 2 , f 3 ), f A = f B = 
f E = max(fi, f 2) f 3 ). In the special case, Jt ->• 0, Eq. ffTSajl 
reduces to a second-degree equation whose solutions are Yoi — Y02 = 1/3, 
while the related function is otherwise divergent as a — > +00. In general, 
the extremum points of VDW isothermal curves ( JT < 1) occur at Y = Yb 
(minimum) and Y = Yd (maximum), Yb < Yd- As ft — > 0, Yb ~ ► 1/3, 
— >■ +00, where, in all cases, l/3<ye<l< Yd- 

The two areas defined by the intersection of a generic VDW isothermal 
curve (JT < 1) and related real isothermal curves (see Fig.[3]), are expressed 
as: 



r vc rvc rye 

W x = I PcdV- pdV = p c V c MTc-Ta)- / #df 
jy A jva JYa 

W 2 = pdV- Pc dV = p c Vc / ^y-MTe-Yc) 

J Vr J Vc. J Vc 



(24a) 
(24b) 



■Fe 

'v c - jv c ~ y-ifc 

and the substitution of Eq. (1141) into ()24|) allows explicit expressions for the 
integrals. The result is: 

W 1 



PcV c 



c(fc - fA) '3 yin 3yw + f A fc ; (25a) 
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J* = | riB |^-l_3(^-ra . (25b) 

PoVfc 3 3f c -l fcfE 

and the condition, W\ = W 2 , after some algebra reads: 
j 8 T , 3f E -l 3 

fc = iywi k iyw"Mi 1 (26) 

where, for a selected isothermal curve, the unknowns are i>c = i>A = i>E, Ta, 
and Ye- 

The reduced volumes, Ya, Yc, Ye, see Fig-El ma Y be considered as inter- 
sections between a VDW isothermal curve ( JT < 1) and a horisontal straight 
line, jf) = $c, in the (OYlft) plane. In other words, J^a, Yc, Ye, are the real 
solutions of the third-degree equation: 

r -^ + \£) r + Tc r -Tc^-' (27) 

which has been deduced from Eq. (1141) . particularized to $ = tfc- The re- 
lated solutions may be calculated using Eqs. (1231 . The last unknown, tfc, is 
determined from Eq. (1261) . 

An inspection of Fig.[3] shows that the points, A and E, are located on the 
left of the minimum, B, and on the right of the maximum, D, respectively. 
Keeping in mind the above results, the following inequality holds: Ya < 
Yb < 1 < Yd — Ye, which implies further investigation on the special case, 
Yc = 1- The particularization of the VDW equation of state, Eq. (fl4~l) . to the 
point, C = Ci, assuming Yc x — 1> yields: 

r = h 4 1 -, (28) 



and Eq. (127|) reduces to: 



f 3 -(l + 26)f 2 + 36f -6 = ; (29a) 

h = ^r ; ( 29b ) 

with regard to the generic third-degree equation, Eq. (fl9~j) . the three solutions, 
Xi, X2, X3, satisfy the relations (e.g., Spiegel, 1968, Chap. 9): 

%i + x 2 + £3 = — oi ; (30a) 

xix 2 + x 2 x 3 + X3X1 = a 2 ; (30b) 

X1X2X3 = -a 3 ; (30c) 
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where, in the case under discussion: 



a\ = — 1 — 2b ; = 36 ; = — b ; (31a) 

xi = f A ; x 2 = f Cl = 1 ; x 3 = f E ; (3lb) 

and the substitution of Eqs. (13~TI) into two among (13"U|) yields: 



f A = 6-v^2T5 ; (32a) 
f E = 6 + v / ^ r 6 ; (32b) 

and the combination of Eqs. ( 128|) . ( 129bj) . and (132]) produces: 



1 -2J1 - T 

n= 4 y_ 3 ; ; (33a) 



1 + 2J1- T 

ft= 4 y_ 3 ; r<i ; ; (33b) 

which, together with = 1, are the abscissae of the intersection points be- 
tween a selected VDW isothermal curve in the {OY'f) plane and the straight 
line, jf) = $c\, in the special case under discussion. 

The substitution of Eqs. (133]) into (|26|) . the last being related to the real 
isothermal curve, yields: 



In V . =6 ; (34) 



1- f 3-2^-3^1- f 

which holds only for the critical isothermal curve, = 1. Accordingly, the 
abscissa of the intersection point, C, between a selected VDW isothermal 
curve and related real isothermal curve, see Fig.[3l cannot occur at fc — 1 
unless the critical isothermal curve is considered. Then the third-degree equa- 
tion, Eq. f[2T|) . must be solved in the general case by use of Eqs. (123]) . The 
results are shown in Tab.[H where the following parameters (in reduced vari- 
ables) are listed for each VDW isothermal curve, see Fig. El the temperature, 
JT, the lower volume limit, Ya, for which the liquid and vapour phase coexist; 
the extremum point (minimum) volume, fs', the intermediate volume, ^c, 
for which the pressure equals its counterpart related to the corresponding 
lower and upper volume limit, for which the liquid and vapour phase coexist; 
the extremum point (maximum) volume, )^b; the upper volume limit, f^, for 
which the liquid and vapour phase coexist; the extremum point (minimum) 
pressure, j/>b] the pressure, j/>a = i>c = i>Ei related to the horisontal real 
isothermal curve; the extremum point (maximum) pressure, $d. The locus 
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Table 1: Values of parameters, f, f A , f B , f c , Yd, Ye, i>B, ic, i>D, within 
the range, 0.85 < < 0.99, using a step, A ft = 0.01. All values equal unity 
at the critical point. Index captions: A, C, E - intersections between VDW 
and real isothermal curves; B - extremum point of minimum; D - extremum 
point of maximum. Extremum points are related to VDW isothermal curves, 
while their real counterparts are flat in presence of both liquid and vapour 
phase. For aesthetical reasons, 01 on head columns stands for unity. 
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8.4584 
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8 


0593 
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0336 


1.2867 


1.6118 


7.9811 


8.4762 


8.7319 


0.97 


7.3756 


8 


2830 


1 


0249 


1.2404 


1.4960 


8.5328 


8.8429 


9.0185 


0.98 


7.7554 
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5611 
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0164 


1.1892 


1.3761 


9.0576 


9.2191 


9.3209 


0.99 


8.3091 
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9461 


1 


0081 


1.1278 


1.2430 


9.5510 


9.6048 


9.6437 
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of the intersections between VDW and real isothermal curves is represented 
in Fig. [2] as a trifid curve, where the left, the right, and the middle branch 
correspond to ^a, and fc, respectively. The common starting point 
coincides with the critical point. The locus of the VDW isothermal curve 
extremum points is represented in Fig.|2]as a dotted curve starting from the 
critical point, where the left and the right branch corresponds to minimum 
and maximum points, respectively. 

A fluid state can be represented in reduced variables as (f, j/>, JT), where 
one variable may be expressed as a function of the remaining two, by use 
of the reduced ideal gas equation of state, Eq. (fl2"]h or the reduced VDW 
equation of state, Eq. ffl3|) . The formulation in terms of reduced variables, 
Eqs. (11 ip . makes the related equation of state universal i.e. it holds for any 
fluid. Similarly, the Lane-Emden equation expressed in polytropic (dimen- 
sionless) variables, describes the whole class of polytropic gas spheres with 
assigned polytropic index, in hydrostatic equilibrium (e.g., Chandrasekhar 
1939, Chap. IV, §4). 

The states of two fluids with equal ( f ■, i>i T)i are defined as corresponding 
states. The mere existence of an equation of state yields the following result. 

Law of corresponding states. Given two fluids, the equality between two 
among three reduced variables, f , jf>, fl, implies the equality between the 
remaining related reduced variables i.e. the two fluids are in corresponding 
states. 

The law was first formulated by van der Waals in 1880. For further details 
refer to specific textbooks (e.g., LL67, Chap. VIII, §85). 



3 Astrophysical fluids 

Let macrogases be defined as two-component fluids which interact only grav- 
itationally. The virial theorem for subsystems reads (Caimmi et al., 1984; 
Caimmi and Secco, 1992; CV08): 

2(E u ) kiQ + (E uv ) VlT = ; u = i,j ; v = j,i ; (35a) 
{E uv ) v i r = (E u ) se \ + (E uv ) t i<i ; (35b) 

where i and j denote the inner and outer subsystem, respectively, E kin is 
the kinetic energy, E seh E tid , and E vir , are the self, tidal, and virial potential 
energy, respectively. The related definitions are: 

If 3 

(E u ) kin = - pu(x 1 ,x 2 ,x 3 )J2( v u) 2 s d 3 S u ; (36) 
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(E u ) sc i= I pu(xi,X2,X 3 )J2 X s^-^ 3 Su 
J S u S=l 5 

1 r „ 
= ~o Pu(x 1 ,x 2 ,x 3 )V u (x 1 ,x 2 ,x 3 ) d S u ; (37) 
^ J s u 

f 3 dV v 

(E uv ) tid = Pu(x 1 ,x 2 ,x 3 )J2 x s-K-^d 3 S u ; (38) 

where p is the density, i> s the velocity component, S the volume, and V the 
gravitational potential. 

The virial theorem makes a necessary (but not sufficient) condition for 
dynamical or hydrostatic equilibrium, which implies the parameters of the 
virialized configuration must be considered as averaged on a sufficiently long 
time. For further details refer to specific textbooks (e.g., Landau and Lif- 
chitz, 1966, Chap. II, §10) and to an earlier paper (Caimmi, 2007). On 
the other hand, general trends exhibited by virialized configurations hold, 
in particular, for self-consistent density profiles implying nonnegative dis- 
tribution functions. Accordingly, density profiles shall be chosen regardless 
from their self-consistency, aiming to investigate general trends instead of 
local properties. To avoid the determination of the gravitational potential, 
which is the most difficult step towards an explicit formulation of poten- 
tial energies, future considerations shall be restricted to homeoidally striated 
ellipsoids (Roberts, 1962). The following results are taken from earlier at- 
tempts (CV08, and further references therein), to which an interested reader 
is addressed. 

The isopycnic (i.e. constant density) surfaces are defined by the following 
law: 

pu = plfu{£u) ; /u(l) = 1 ; u = i,j ; (39a) 
<vu = — p j < £ u < c u ; ^ u = — ; (39b) 

where the dagger denotes a selected reference isopycnic surface, r is the radial 
coordinate along a selected direction, R the related truncation radius, £ a 
related scaled radial coordinate, and S the related scaled truncation radius. 
The mass and the self potential energy read: 

M u = K) mas Mt ; Ml = ^pl(a u )\(a u )i(a u )l ; (40a) 

K) mas = - / F u ^ u )d^ u ; F U (Q = 2 / UUZu d£ u ; (40b) 
2 Jo J(u 

G(M r ) 2 9 /" s « o 

(E u ) scl = -(v u ) scl -^-B u ; (u u U = — F^ u )d£ u ; (41) 
(al)i 16 Jo 



15 



where a\ are semiaxes of the reference isopycnic surface, u mas and u se i are 
profile factors which depend only on the mass distribution, and B is a shape 
factor. For homogeneous configurations, z/ mas = S 3 and u sc \ = (3/10) S 5 . For 
spherical shapes, B = 2. 

Under the further restriction of similar and similarly placed boundaries, 
the following relations hold: 



Z/ T (^i)mas ™ T 



W = ; 2/ - 4 5 m = -^ ; 771 = — 4 : (42b) 
which makes tidal and virial potential energy reduce to: 

(.tiuv) xxx 7 +t (^tra)xxx-^ j (43a) 

u — i,j; v — j,i ; xxx = tid, vir ; (43b) 
and the explicit expression of the profile factors reads: 

(^id = -^mV ext )(r?) ; (vfiU = -l^w<M(rj) ; (44a) 

O O 777,' 

(^w)vir = (^«) S ei + (j / uv )tid ; u = i,j ; v=j,i ; (44b) 
r/ = ^ = ^ ; y>l ; (44c) 

where the functions, iy( mt ) and t// ext \ are defined as: 



M) / 



" dF, 



(77) = ^ F,(0)^0d^ ; (45a) 

^ (ext) (^) = jf ^te)^0d0 ; (45b) 

for further details refer to Appendix |A] In conclusion, Eqs. (l3~9]) - (j4"5T) allow the 
calculation of the virial potential energy for homeoidally striated ellipsoids 
related to similar and similarly placed boundaries. 
The fractional virial potential energy, reads: 



(Eji) viT _ (m f ) 2 (Uji) vir _ 777? 
(iSy)vir V ] (^ij)vir y 



-\ 2 



i ) mas 



; (46) 



which, for assigned density profiles, depends on either the reference fractional 
mass, 777', and the fractional scaling radius, y\ or the fractional mass, m, 
and the fractional truncation radius, y. 
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Strictly speaking, Eq. ( 14"6|) is valid provided the indices, i and j, denote the 
embedded and the embedding subsystem, respectively, which implies y > 1. 
If the role of the two subsystems is reversed, < y < 1, it has to be kept 
in mind that the inner and the outer component are denoted by the indices, 
j and i, respectively. Then the quantities of interest must be calculated 



{m)) l ,y y 1 , 



y ] ->■ (y ] 



according to the changes, m — > mT 1 . 
i -H- j, which yields <fi = (i%) v ir/(-%i)vir in the domain, y > 1. Following the 
above mentioned procedure where, in addition, <fi — > <f> , allows the explicit 
expression of the fractional virial energy, (f> = {Eji) v - ir / (i^) v i r , in the domain, 
< y < 1, which extends the whole domain to < y < +oo. 

In absence of truncation radius, 5 — > +oo, r\ — > +oo, the reversion occurs 
when the density drops to zero and nothing changes except in infinitesimal 
terms of higher order and infinite terms of lower order. Accordingly, there is 
no need to perform the reversion in this case. 

The combination of Eqs. (1441) and (TJ6]) yields: 



t)^(^)sel-^" (mt) ^ 



j /mas 



2 o \ V i)se\ 
Tfl 



9 5j (l/jQmas y (int) 



9 ZA 



(47a) 



(47b) 



i/mas (ext) 

-mw ( 1 



(v) 



j /mas 



where, for assigned scaled truncation radii, Hj and Ej, the independent vari- 
ables are m^, yt, or m, y. In terms of a second-degree equation in or m, 
Eqs. flUD read: 



Ax 2 + Bx + C = 



9 

B = — k B 
8 

C = -/c C (^)sel 



A; A = 1, ^ 



&c = y\ y ; 



] / mas 



(48a) 
(48b) 

(48c) 
(48d) 



where the positive solution is: 



x 



B 2 - AAC - B 
2A 



x 



(49) 
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and the negative solution has been disregarded due to the lack of physical 
meaning. 

In the special case of coinciding density profiles, fi = fj, Fi = Fj, and 
scaled truncation radii, Sj = S 3 -, fractional masses and truncation radii also 
coincide, m) = m, y^ = y, via Eqs. (142ft . and the same holds for the pro- 
file factors, 0i) mas = (fj)mafl, (^)sci = O j)sei, which depend on the scaled 
truncation radii. Accordingly, Eqs. (j47aj) . ( 14 7b p . also coincide. 

To get a closer analogy with real gases, let Eqs. (jl7|) be rewritten after a 
change of variables, as: 

9 1 w( int \r]) 



x p x v ^fck^W = 



(50a) 



8 (^i)sel 

9 Site)™* 1 u/^fa) 



O (^ijsel 

r; = ~,X V = H.-Xy ; y > 1 ; (50c) 
Xt = (mt) 2 ; Xt=l ; Xt=0 ; tft( Hi , Sj -) = ^ ; (50d) 



v J mas 



(50e) 



(fj)sel 



X p = m 2 ; X v = - ; X T = ; X(H;, S,-) = ^ 

which can be conceived as an equation of state for macrogases. The vari- 
ables, Xy, X p , X T , play a similar role as the volume, the pressure, and the 
temperature, for ordinary gases. Accordingly, X v , X p , X T , shall be defined 
as macrovolume, macropressure, and macrotemperature, respectively. 
The combination of Eqs. (T42|) . ( I50dj) . and fl50e[) yields: 

-i 2 



Xy — — Xy ; Xl 



p 



X p ; X T — X T ; (51) 

V J / mas _ 

which links the variables, (X v , X|,X^), to (Xv,X p ,Xt), and vice versa. 

Strictly speaking, the macrogas equation of state should be deduced from 
dimensional (instead of dimensionless) virial equations for an assigned sub- 
system, as outlined in Appendix [B] On the other hand, a description in 
terms of dimensionless variables turns out to be more useful. 

If the interaction terms are omitted, w/ int ) = u>( ext ) = 0, Eqs. (I50al) and 
(I50bj) reduce to: 

*t*t =Ar t( SijSi )j4 . (52a) 
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(52b) 



which may be considered as equation of state of ideal macrogases, where 
"ideal" means "the interaction terms are omitted". 

The parameters, K* and K, appearing in either macrogas equation of 
state, depend on the scaled truncation radii, Ej and Sj, and on the selected 
density profiles. In other words, the macrogas equation of state is not uni- 
versal, but takes a different form for different density profiles. A restricted 
number of special cases shall be studied below, grounding on earlier results 
(CV08), to which an interested reader is addressed for further details. In 
any case, the following method shall be used: (i) select two density profiles; 
(ii) fix related scaled truncation radii, Sj and 5^; (iii) choose a macrotem- 
perature, 0; (iv) plot related macroisothermal curves, by solving Eq. (150aj) or 

(TO. 

3.1 UU macrogases 

The related density profiles maintain uniform, which is equivalent to poly- 
tropes with index, n = (e.g., Chandrasekhar, 1939, Chap. IV, §4; Caimmi, 
1986), but implies negative distribution functions for stellar fluids (Vander- 
voort, 1980). The particularization of the general expressions to the case 
under discussion, yields for the quantities of interest (CV08): 



U{€u) = 1 ; < £ u < E u ; u = i,j ; (53) 

F u (£ u ) = ~l - £ ; u = i,j ; (54) 

K)mas = 23 . U = iJ ; (55) 

3 

K)sei = — ^1 ; u = i,j ; (56) 

^(Jl) = ~^^-l) ; (57) 

w^\ V ) = ■ (58) 

where Eqs. (157]) and fl58l hold under the additional restriction of similar and 
similarly placed boundaries (CV08). 

In the case under consideration of uniform density profiles, without loss 
of generality, it can be assumed a scaled truncation radius, E u = R u /rl = 1, 



which implies y = y\ m = m\ due to Eqs. (T42i) and (155]) . Accordingly, 
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Eq. (US]) reduces to (CV08): 



(tttJ) 2 / y 



5 1 + 



W 1 ^5 



2* 



1 + 



m 1 



2/>l 



(59a) 



77V (t/) 2 



1 + 



777' 



,t !. 



0<y<l ; 



2 2^ 



(59b) 



where y is the outer to inner ellipsoid axis ratio, Rj/ Ri, according to Eq. (142b)) . 
In addition, Eqs. (|48l) and (H91 reduce to: 



m = \J f3 2 + y4> — ft 



1 i 

2 7/2 

A = — 
10 



3 
2 



c = — # 
io yr 



y > i 



< 7/ < 1 

5 



(60a) 
(60b) 

(60c) 

(60d) 



where the negative solution is not considered due to the lack of physical 
meaning. 

The explicit expression of the square fractional mass, m 2 , extracted from 



Eq. (I60ap in dimensionless variables, X p = = m 2 
l/ y = l/ y \ X T = X ] T = 0, Eqs. fl50dj) . f lBDejl . reads: 

X„ 



777 



t) 2 , X 



v 



X- 



X, 



2/3 2 + -i- - 2/3W/3 2 + -i- 



X 



2/3 = X 



v 

2 /5 1 3 
v UxJ"" 2 



v 



X, 



or, more explicitly: 



X n 



1,2 X 2 



-X 



v 



5 1 3 
2X$~ 2 



-X v - - ) X T 



< X v < 1 
X v > 1 



(61a) 
(61b) 
(61c) 



V 2X 2 



Xn 



Xn 



+ 



Xj 
X 



v 



1/2 



< X v < 1 



(62a) 
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which is the actual UU (AUU) macrogas equation of state. 

The ideal UU (IUU) macrogas equation of state is obtained by the com- 
bination of Eqs. PUD , (j52bj) , fl55]), and (JH]). The result is: 



*p = 1T ! (63) 
Ay 

which represents a hyperbola with equal axes, for fixed Xt- 

Macroisothermal curves related to IUU (tidal potential energy excluded) 
and AUU (tidal potential energy included) macrogases, are plotted in Fig.lU 
left and right panel, respectively, for values of the macrotemperature, Xt = 
0.85, 0.90, 0.95, 1.00, 1.05, 1.10, from bottom to top. The coordinates, X v = 
X v , X p = XjJ, Xt = X T , may be conceived as normalized to their fictitious 
critical counterparts, X Vc = Xy c = 1, X Pc = X+ = 1, X Tc = = 1, as 
<p = m = m) for y — y* — 1, according to Eqs. (1591) or (11231) . which implies 
= 1 for m = m) = 1. The comparison with ideal and VDW gases, plotted in 
Fig.[lJ shows a similar trend, except the absence of a critical macroisothermal 
curve, above which the extremum points disappear. 

Contrary to ordinary gases, no experiment can be performed on macro- 
gases to ascertain the existence of a phase transition moving along a selected 
macroisothermal curve, where the path is a horisontal line instead of a curve 
including the extremum points. Then the existence of the above mentioned 
phase transition and flat real macroisothermal curves, must necessarily be 
assumed as a working hypothesis, by analogy with VDW isothermal curves 
(below the critical one). The loci of extremum points of AUU macroisother- 
mal curves plotted in Fig.H] (right panel), are represented as dotted lines in 
Fig.EJ 

Unlike the VDW equation of state, Eq. ffT^j) . the AUU macrogas equation 
of state, Eq. (1<32~j) . is not analytically integrable. Then the procedure used 
for determining a selected macroisothermal curve, must be numerically per- 
formed. The main steps are (i) calculate the intersections, Xy A , Xy c , Xy E , 
Xy A < Xy c < Xy E , between the generic horisontal line in the (OXyX p ) 
plane, X p = const, and the AUU macrogas equation of state, within the 
range, X PB < X p < X PD , where B and D denote the extremum points of 
minimum and maximum, respectively; (ii) calculate the area of the regions, 
ABC and CDE; (iii) find the special value, X p = X po , which makes the 
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two areas equal; (iv) trace the real UU (RUU) macroisothermal curve as 
a horisontal line connecting the points, (X Va ,X pa ), (X Vc , X pc ), (X Ve , X PE ), 
Xp A = X pc = X PE = X Po . 

The loci of the intersections between AUU and RUU macroisothermal 
curves, are represented as dashed lines in Fig. [51 The absence of a critical 
macroisothermal curve makes a band-like instead of a bell-shaped region exist 
on the (OXyX p ) plane, as shown by comparison with Fig.[2J The AUU and 
RUU macroisothermal curves are plotted in Fig.[6l with regard to the special 
case, Xt = 0.85. Values of parameters, X?, Xy A , Xy B , Xy c , Xv D , Xy E , X PB , 
X Pc , X PD , are listed in Tab. [2] within the range, 0.85 < Xt < 1.10, using a 
step, AX T = 0.01. 



3.2 HH macrogases 

The related density profiles exhibit a central cusp and null value at infinite 
distances (Hernquist, 1990), and have been proved to be consistent with 
nonnegative distribution functions, in an acceptable parameter range (Ciotti, 
1996). The particularization of the general expressions to the case under 
discussion, yields for the quantities of interest (CV08): 

g 

/u(CJ = c 7z , c no ; < i u < E u ; u = i,j ; (64) 

W = (r^-(W ; u = h3 ; (65) 

12S 2 

u = i,j ; (66) 



(1 + H U ) 2 
12S|(4+ ~ u 
T+H 



fruU = 7, ; u = i,j ; (67) 



w M)( 77 ) = _1282/ t (i- ' 



2(yt 



(yt _ l)Vr? + 2(yt - 1)^ 



+ (yt-i)(yt + % t ?7+2(V + i)in 



(j/fy + I) 2 
?7 + 1 



1 + r/ 



?/l?7 + 1 
1 1 1 



y^rj + 1 



2y^r] + 1 

2 (1 + H,-) 2 (yt)a f _ (y^ + l)3j 



u;< te >fa) = -128 - 



4/7+1 
(?7 + l) 4 



+ 1 



(68a) 



2(l + S,) 2 (^+i: 



1 ; 



(68b) 
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Table 2: Values of parameters, X^, Xy A , Xy B , Xy c , Xy D , Xy E , X PB , X pc , X pv 
(to be conceived as normalized to their fictitious critical counterparts, X Vc = 
1, X Pc = 1, X Tc = 1), within the range, 0.85 < X T < 1.10, using a step, 
AX T = 0.01. Index captions: A, C, E - intersections between AUU and RUU 
macroisothermal curves; B - extremum point of minimum; D - extremum 
point of maximum. Extremum points are related to AUU macroisothermal 
curves, while their RUU counterparts are flat within the range, Xy A < Xy < 
Xy E . For aesthetical reasons, 01 on head columns stands for unity. 
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using Eqs. (jHJ) and fl66l) - (|69|) . the actual HH (AHH) macrogas equation of 
state is obtained from the particularization of Eqs. (15 Op to the case of interest 
for the domain, y > 1. The extension to the domain, < y < 1, can be done 
following the procedure outlined above in dealing with Eq. (I46p . 

In absence of truncation radius, H — > +00, r\ — > +00, and Eqs. f l65|) -( |69|) 
reduce to (CV08): 
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using Eqs. and f!7T]) - fl71|) . the AHH macrogas equation of state in the 
special situation under discussion, is obtained from the particularization of 
Eqs. (150]) to the case of interest. 

The ideal situation, where the interaction terms are omitted, is obtained 
using Eqs. (152]) instead of Eqs. (150]) . More specifically, the ideal HH (IHH) 
macrogas equation of state is derived by the combination of Eqs. (I50ej) . (]52bl) . 
(lolj]) . and (16T]) . The result is: 

which represents a hyperbola with different axes (unless Sj = Hj), for fixed 
At. 

Macroisothermal curves (^T p = X p /X Pc vs. Jfiy = Xy/X Vc ) related to 
IHH (tidal potential energy excluded) and AHH (tidal potential energy in- 
cluded) macrogases, are plotted in Fig. [7] left and right panels, for different 
values of scaled truncation radii, (H^Hj), labelled on each panel, and same 
values of the reduced macrotemperature, JjC^ = X T /X Tc = 0.90, 0.95, 1.00, 
1.05, 1.10, 1.15, from bottom to top. The limit, (Hj,Ej) —> (+oo,+oo), 
makes only little changes. The comparison with ideal and VDW gases, plot- 
ted in Fig. [I] shows a similar but reversed trend. More specifically, extremum 
points occur above, instead of below, the critical macroisothermal curve. A 
complete analogy can be obtained using the transformations, Xy — > 1/Xy, 
X p -> l/X p , X T -> 1/X T . 

The existence of a phase transition moving along a selected macroisother- 
mal curve, where the path is a horisontal line instead of a curve including 
the extremum points, must necessarily be assumed as a working hypothesis, 
due to the analogy between VDW isothermal curves and AHH macroisother- 
mal curves. As in the case of UU macrogases, AHH macroisothermal curves 
must be numerically determined, following the same procedure outlined in 
Subsection 13. 1L Characteristic loci of AHH macroisothermal curves plotted 
in Fig. [7J (right panels), are represented in Fig. [8] The loci of intersections 
between AHH and real HH (RHH) macroisothermal curves, are represented 
in Fig. |H] as trifid curves, where the left branch corresponds to ^Ty A , the right 
branch to ^V E , and the middle branch to fiv c - The critical point is the 
common starting point. The loci of AHH macroisothermal curve extremum 
points are represented in Fig. [8] as dotted curves starting from the critical 
point, where the left branch corresponds to minimum points and the right 
branch to maximum points. The RHH macroisothermal curves, when differ- 
ent from their AHH counterparts, lie within the larger bell-shaped regions. 
The limit, (Sj,Hj) — > (+oo,+oo), makes only little changes. Values of pa- 
rameters, ^Tt, ftv A , $v B , ^V c , ^V D , ftv E , ^ PB , fcp c -, ftp D ' are listed in Tab. [3] 
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within the range, 1.00 < ^Tt < 1.15, using a step, A^T T = 0.01, in the 
limit, (Hf, Hj-) — > (+oo,+oo). The AHH macroisothermal curves are plotted 
in Fig.[9]for scaled truncation radii, (Hi, Hj), as in Fig.[BJ with regard to the 
special case, ^tr = 1.15. The limit, (Hi, Hj) — » (+oo,+oo), makes only little 
chang shown by comparison with the dashed curve. 

Values of parameters, m^, m, y^ , ?/, 0, related to the critical macroisother- 
mal curve, for selected scaled truncation radii, (H^H,), are listed in Tab.HJ 
A better method has been used with respect to an earlier attempt (CV08), 
yielding improved results. 



3.3 HN/NH macrogases 



A description of H density profiles has been provided in Subsection 13.21 The 
remaining N density profiles also exhibit a central cusp and null values at in- 
finite distances (Navarro et al., 1995, 1996, 1997). Mass distributions defined 
by an inner H and outer N density profile were found to be self-consistent, 
in an acceptable parameter range, with regard to the non negativity of the 
distribution function (Lowenstein and White, 1999), using a theorem stated 
in an earlier attempt (Ciotti and Pellegrini, 1992). The particularization of 
the general expressions to N density profiles and HN macrogases, yields for 
the quantities of interest (CV08): 
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Table 3: Values of parameters, ^T T , ^ty A , ]£ Vb , $ Vc , ^ty D , ^Ty E , fipm fei $pd> 
within the range, 1.00 < ^Tt < 1.15, using a step, A^Tt = 0.01, except near 
the critical point, where convergence was not attained. All values equal unity 
at the critical point. Results are related to infinitely extended configurations, 
(Sj,Sj) — > (+oo,+oo). Index captions: A, C, E - intersections between 
AHH and RHH macroisothermal curves; B - extremum point of minimum; 
D - extremum point of maximum. Extremum points are related to AHH 
macroisothermal curves, while their RHH counterparts are flat within the 
range, ^Ty A < jCy < fiv E - For aesthetical reasons, 01 on head columns 
stands for unity. 
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Table 4: Values of the scaling fractional mass, m\ the fractional mass, m, 
the scaling fractional radius, y\ the fractional truncation radius, y, and the 
fractional energy, 0, related to the critical point i.e. the horizontal inflex- 
ion point on the critical macroisothermal curve, for selected scaled trunca- 
tion radii, (Si, Sj), of HH density profiles. In absence of truncation radius, 
(Si, Ej) — > (+00, +00), the results are independent of y/y^ = Sj/S». 
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using Eqs. fj44"j) . fl66l) . (|67|) . and (T78|) - (l8T|) . the HN macrogas equation of state 
is obtained from the particularization of Eqs. (|50|) to the case of interest for 
the domain, y > 1. The extension to the domain, < y < 1, can be done 
following the procedure outlined above in dealing with Eq. ( 14"6"|) . 

To this aim, the counterparts of Eqs. f lHOj) and ( IHTj) . related to NH macro- 
gases (i = N, j = H), are needed. The particularization of Eqs. fj4"5]) to the 
case under discussion, yields: 
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where Sj = S^, S,- = S#, while the contrary holds with regard to Eqs. ([80]) 
and flH}. Using Eqs. ®, flM}, (EZD , JSj), (JZ2D, QS2), and (JHJ), the NH 
macrogas equation of state is obtained from the particularization of Eqs. ( [50]) 
to the case of interest for the domain, y > 1, which corresponds to < y < 1 
for HN macrogases and vice versa. 

In absence of truncation radius, S — > +00, 77 — » +00, and Eqs. ([77]) - (jST]) 
reduce to: 
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where the self potential-energy profile factor remains finite, although the 
mass profile factor undergoes a logarithmic divergence. Similarly, Eqs. 
and (l83l) reduce to: 
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where Sj = Ejv, £j = while the contrary holds with regard to Eqs. 
and (JEED. Using Eqs. (]4l]h (15B]h (157]h (f7T]), (f72]h and (f78>(j9"0]h the actual HN 
(AHN) and NH (ANH) macrogas equation of state in the special situation 
under discussion, is obtained from the particularization of Eqs. (1501 to the 
case of interest. 

The ideal situation, where the interaction terms are omitted, is obtained 
using Eqs. (152"]) instead of (1501 . More specifically, the ideal HN (IHN) and 
NH (INH) macrogas equation of state is obtained from the combination of 
Eqs. (l50el . fl52b|) . ([66]), (EZJ), (J7HD, and ([79]). The result is: 



X 1 4 + Eh [(l + ^lnq + H.y)-^] 2 X T 
p 3 E H E N (2 + E N )-2(l + E N )\n(l + E N )X v ' 1 J 

which represents a hyperbola with different axes for fixed Xt- In the limit of 
infinite extension, E — >• +oo, both the left and right-hand side of Eq. (19T]) di- 
verge, but a different equation of state may be derived starting from Eqs. ( I50dj) 
and (|52ap . following a similar procedure. The result is: 



x t = 1 5^(4 + 5 g ) {l + E N f X} 

p 3 {1 + EhY E N (2 + E N )-2(l + E N )\n(l + E N )x\ J 



(92) 



which also represents a hyperbola with different axes for fixed X T . 

Macroisothermal curves (X p = X p /X Pc vs. Xy = X v /XvJ related to 
IHN/INH (tidal potential energy excluded) and AHN/ANH (tidal potential 
energy included) macrogases, are plotted in Fig. [TO] left and right panels, for 
different values of scaled truncation radii, (Hi, Ej), labelled on each panel, and 
same values of the reduced macrotemperature, Xt = Xt/Xt c = 0.90, 0.95, 
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1.00, 1.05, 1.10, 1.15, from bottom to top. The limit, Ej) — > (+oo, +oo), 
makes only little changes. The comparison with ideal and VDW gases, plot- 
ted in Fig.[lJ shows a similar but reversed trend. More specifically, extremum 
points occur above, instead of below, the critical macroisothermal curve. A 
complete analogy can be obtained using the transformations, Xy — > 1/Xy, 
X p -> 1/X P , X T 1/X T . 

The existence of a phase transition moving along a selected macroisother- 
mal curve, where the path is a horisontal line instead of a curve including the 
extremum points, must necessarily be assumed as a working hypothesis, due 
to analogy between VDW isothermal curves and AHN/ANH macroisother- 
mal curves. As in the case of UU macrogases, AHN/ANH macroisothermal 
curves must be numerically determined, following the same procedure out- 
lined in Subsection 13.11 The loci of AHN/ANH macroisothermal curves plot- 
ted in Fig.[T0] (right panels), are represented in Fig.[TTJ The loci of intersec- 
tions between AHN/ANH and real HN/NH (RHN/RNH) macroisothermal 
curves, are represented in Fig.[TT]as trifid curves, where the left branch cor- 
responds to ^V A , the right branch to ^ty E , and the middle branch to fiv c - 
The critical point is the common starting point. The loci of AHN/ANH 
macroisothermal curve extremum points are represented in Fig. [11] as dotted 
curves starting from the critical point, where the left branch corresponds to 
minimum points and the right branch to maximum points. The RHN/RNH 
macroisothermal curves, when different from their AHN/ANH counterparts, 
lie within the larger bell-shaped regions. The limit, — > (+oo,+oo), 
makes only little changes. Values of parameters, ^Tt, ftv A , $v B , ^V c > fiv B , 
^Tve, ^T Pb , $p c , fi PD , are listed in Tab. [5] within the range, 1.00 < ^Tt < 
1.15, using a step, A^Tt = 0.01, in the limit, — > (+oo,+oo). The 
AHN/ANH macroisothermal curves are plotted in Fig. [15] for scaled trunca- 
tion radii, (S^EL,-), as in Fig.[TTJ with regard to the special case, ^Tt = 1.15. 
The limit, (Sj,Hj) — > (+oo, +oo), makes only little changes, as shown by 
comparison with the dashed curve. 

Values of parameters, m) , m, y\ y, <ft, related to the critical macroisother- 
mal curve, for selected scaled truncation radii, (Sj, Sy), are listed in Tab. [61 
A better method has been used with respect to an earlier attempt (CV08), 
yielding improved results with the occurrence of the critical macroisothermal 
curve in all cases. In absence of truncation radius, (Sj,Sj) — > (+oo, +oo), 
small differences appear in the results, probably due to (i) the divergence of 
the fractional mass, m — > +oo, as H^r — > +oo, and (ii) numerical calculations 
have been performed, taking H ~ 10 10 . Related results listed in Table [6] may 
be considered as weighted means. 
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Table 5: Values of parameters, ^Tt, ftv A , ftv B , ftv c -> ftv u , ftv E , ft PB , fi Pc , 
ft PD , within the range, 1.01 < ^Tt < 1.15, using a step, A^Tt = 0.01. All 
values equal unity at the critical point. Results are related to infinitely ex- 
tended configurations, (Sj,Sj) — > (+oo, +oo). Index captions: A, C, E - 
intersections between AHN/ANH and RHN/RNH macroisothermal curves; 
B - extremum point of minimum; D - extremum point of maximum. Ex- 
tremum points are related to AHN/ANH macroisothermal curves, while their 
RHN/RNH counterparts are flat within the range, ^Ty A < J/Cy < ^ty E . For 
aesthetical reasons, 01 on head columns stands for unity. 
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Table 6: Values of the scaling fractional mass, nv, the fractional mass, m, 
the scaling fractional radius, y\ the fractional truncation radius, y, and the 
fractional energy, 0, related to the critical point i.e. the horizontal inflexion 
point on the critical macroisothermal curve, for selected scaled truncation 
radii, (Sj,Sj), of HN/NH density profiles. In absence of truncation radius, 
(Sj, Ej) — > (+00, +00), the results are independent of y/y^ = Sj/Sj. 
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4 Application to elliptical galaxies 

The luminosity-weighted second moment of the line-of-sight velocity distri- 
bution within the half-light radius has recently been determined for samples 
of early-type galaxies, using integral-field spectroscopy such as SAURON 
(SIV; SX). Compared to the central velocity dispersion, which was some- 
times used before, the above mentioned quantity has the advantage that it 
is only weakly dependent on the details of the aperture used. In addition, 
it is an approximation to the second velocity moment which appears in the 
virial equations (e.g., Binney and Tremaine, 1987, Chap. 4, §4.3). In other 
words, it may be conceived as a rms mass-weighted velocity, with a weak 
dependence on the features of the orbital distribution. For further details 
refer to the parent papers (SIV; SX). 

An application of the current model to SAURON sample objects can be 
performed along the following steps: (i) select SAURON data of interest; (ii) 
calculate parameters appearing in the virial equations; (iii) make a corre- 
spondence between model galaxies and sample obiects; (iv) represent model 
galaxies as points on the (O^Ty^p) plane. 
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4.1 Data selection 



The sample used (CV08, N = 16) is made of elliptical galaxies, extracted 
from richer samples of early-type galaxies investigated within the SAURON 
project (SIV, N = 25; SX, N = 48). The selection has been restricted 
to elliptical galaxies common to the above mentioned samples for two main 
reasons, namely (i) the current model best applies to elliptical galaxies and 
their hosting haloes, and (ii) some physical parameters of interest are listed 
in either S IV or S X. The whole data set needed for the application is shown 
in Tab. [3 For further details refer to the parent papers (SIV, columns 2-6; 
SX, columns 7-11). 

Values to be actually used in the application are related to: the effective 
(half-light) radius, R e (S IV); the total observed I band galaxy magnitude, It 
(SIV); the mass- luminosity ratio of the stellar population, Mi/L (SIV); the 
galaxy distance modulus, m — M (hats avoid confusion with the fractional 
mass, m = Mj/Mj, and the total mass, M = Mi + Mj) (S IV); the luminosity- 
weighted average ellipticity, ei, on a plane perpendicular to the line of sight, 
within either an isophote enclosing an area, A = nR^, or the largest isophote 
fully contained within the SAURON field, whichever is smaller (SX); the 
luminosity- weighted squared mean velocity component, parallel to the line 
of sight, within either an ellipse of area, A, ellipticity, e±, and related posi- 
tion angle, or the largest similar ellipse fully contained within the SAURON 
field, whichever is smaller (SX); the luminosity- weighted squared velocity 
dispersion, parallel to the line of sight, within either an ellipse of area, A, 
ellipticity, e±, and related position angle, or the largest similar ellipse fully 
contained within the SAURON field, whichever is smaller (SX). 

Parameters listed to gain further insight even if not used in the application 
are: the ratio between maximum radius, R max , sampled by the kinematical 
observations, and effective radius, R e (SIV); the ratio between square root 
luminosity-weighted squared mean velocity component and dispersion veloc- 
ity component, parallel to the line of sight (SX); the kinematic classification 
(F - fast rotator; S - slow rotator) (SX). 

For further details refer to the parent papers (S IV; S X) and an earlier 
attempt (Binney, 2005). 

4.2 Determination of model parameters 

With regard to the stellar subsystem, two model parameters can directly 
be inferred from the data (CV08). More specifically, stellar masses are de- 
duced from luminosities and mass-luminosity ratios (in / band), as Mj/Mio = 
(L/L )[(M,/L)/(lO lo m /L )]; L/L = exp 10 {-0.4[/ T — (m — M) — 4.11]}; 
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Table 7: Data related to a sample (N = 16) of elliptical galaxies, extracted 
from larger samples of early-type galaxies investigated within the SAURON 
project (SIV, N = 25; SX, iV = 48), which are used in the current pa- 
per. Column captions: (1) NGC number; (2) effective (half-light) radius, R e , 
measured in the I band (SIV); (3) ratio between maximum radius, R max , 
and effective radius, R e (SIV); (4) total observed I band galaxy magnitude 
(SIV); (5) mass-luminosity ratio of the stellar population (SIV); (6) galaxy 
distance modulus (hats avoid confusion with the fractional mass, m, and the 
total mass, M) (SIV); (7) luminosity-weighted average ellipticity, e±, on a 
plane perpendicular to the line of sight, within either an isophote enclos- 
ing an area, A = irR^, or the largest isophote fully contained within the 
SAURON field, whichever is smaller (SX); (8) luminosity- weighted squared 
mean velocity component, parallel to the line of sight, within either an ellipse 
of area, A, ellipticity, e±, and related position angle, or the largest similar 
ellipse fully contained within the SAURON field, whichever is smaller (SX); 
(9) luminosity-weighted squared velocity dispersion, parallel to the line of 
sight, within either an ellipse of area, A, ellipticity, e±, and related position 
angle, or the largest similar ellipse fully contained within the SAURON field, 
whichever is smaller (SX); (10) ratio between the square root luminosity- 
weighted squared mean velocity component and dispersion velocity compo- 
nent, parallel to the line of sight (SX); (11) kinematic classification (F - fast 
rotator; S - slow rotator) (SX). For further details refer to the parent papers 
(SIV; SX) and an earlier attempt (Binney, 2005). 
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and scaling radii are calculated from effective radii (in arcsec) and distances, 
asrj/kpc = (i? e /kpc)/f .81 ; i? e /kpc = [(^/arcsec) (d/Mpc)]/206.265; d/Mpc 
= exp 10 [(m — M)/5 — 5]; where the factor, 1.81, is related to an assumed H 
density profile for the inner subsystem, and the factor, 206.265, is related to 
the choice of measure units (CV08). For further details refer to the parent 
paper (SIV). 

Two additional parameters can be inferred by fitting the data with dy- 
namical models. More specifically, the inclination angle, i, is deduced from 
the best fitting two-integral Jeans model (SIV), and the anisotropy parame- 
ter, 5, is determined from the solution of the dynamical models, supposed to 
be axisymmetric (S X). In fact, fast rotators show evidence of large anisotropy 
and axial symmetry, while slow rotators appear to be nearly isotropic and 
moderately triaxial. For further details refer to the parent paper (SX). 

The intrinsic axis ratio, e, is deduced from the computed inclination, 
under the assumption of axisymmetric configurations (S X), using the relation 
(Binney and Tremaine, 1987, Chap. 4, §4.3): 

1 - e 2 = 1 - f hs ; (93) 
sin i 

where e Q b s is the observed axis ratio related to an inclination angle, i, between 
the symmetry axis and the line of sight (i = 90° for edge-on configurations). 

The mass-weighted mean square velocity component and dispersion ve- 
locity component, parallel to the line of sight, for a galaxy observed (obs) at 
an inclination angle, i, under the assumption of axisymmetric (a\ = 02) and 
isotropic on the equatorial plane (an = 022) configurations (SX), are related 
to their edge-on (edo) counterparts as (Binney and Tremaine, 1987, Chap. 4, 
§4.3): 

[< v\ > 1/2 ] obs = [< v\ > 1 / 2 ] cdo sin* ; (94) 
[< v\ > 1/2 ]obs = [< v\ > 1/2 ]cdo(l - 5 cos 2 if 2 ■ (95) 

and the intrinsic mean rotational velocity and velocity dispersion are ex- 
pressed as: 

< >^= V2[< v\ = 4^[< v\ >^} ohs ■ (96) 

< a > ' — (< <T n > + < a 22 > + < (T33 >) 



{2[< a 2 >] cdo + (!-£)[< a 2 >] cdo y/ 2 



1/2 



= (3 - 5f 2 [< a 2 > 1 / 2 ] cdo = ( 3 5 — ) [< a 2 > 1 / 2 ] obs ; (97) 
" \ 1 — cos z 1 1 " 
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in terms of observed quantities, where 5 = 1— < crf 3 > / < o\ x >= 1 — 
< cr| 3 > / < cr| 2 > by definition (e.g., Binney, 2005). 



4.3 Model galaxies vs. sample objects 

The kinetic energy of the stellar subsystem is: 

(E^ n = l -M l {<(^) 2 l > + <a 2 l >) ; (98) 
and the combination of Eqs. (|55]h (HO)), (j3T]h and (19"8]) yields: 

M, {< (a5 « > + < >} = + ^^ffi ; 

1 J [(^)mas] 2 (o})i [(^)mas] 2 (a})l 

(99) 

which, after some algebra, takes the form: 

< <*i > Oi)l 1 f < W)i > + j \ = (^)sel + (^j)tid . (100) 



GM { B { <af> J [(i/, 



2 

« /mas J 



where, for an inner H density profile, the scaling radius, r\ = (aj)i, may be 
chosen as: 

r \ = Yk ; M M) = \m, ; M t {R e )= l -M t ; (101) 

and the shape factor, B, reads (e.g., Chandrasekhar, 1969, Chap. 3 §§17, 22; 
Caimmi, 2009): 



.arcsinvT^ arcsm 



1 - (1- < e_L >) 



B = 2 l^= — = 2 smi = (102) 

/i_(i_<e, >)2 



sin? 

for axisymmetric configurations, where the last equality is owing to Eq. (J93 
and the definition of ellipticity, e = 1 — e. 

The combination of Eqs. dHSD, dHZD, 0, ffTUUD . and (ITUTjl yields: 

[< cr\ >] obs jg e 1 ( 2 [<_gj >]obs 3-j 
2GM i (R e ) 1.81 £ [sin 2 * [< a 2 >] obs + 1 - 6 cos 2 % ' ' " ' ' 

Cij = ^ ± (U f id ; (103b) 
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where the left-hand side of Eq. ( I103al) is expressed in terms of quantities 
which are deduced from either observations or fitting with dynamic models, 
and then may be determined for an assigned sample object. Conversely, 
the right-hand side of Eq. fll03aj) depends on the selected density profiles, 
and then may be determined for an assigned model galaxy. In this view, 
Eq. fll03ap may be read as a correspondence between sample objects (left) 
and model galaxies (right). 

The dimensionless energy, Cy, defined by Eq. fll03b|) . depends on four 



variables via Eqs. (140]) - fl45]) : the scaled truncation radii, Sj, Ej, the fractional 
mass, rrv (or m), and the fractional radius, y^ (or y). 

4.4 Model galaxies on the (O^ty^p) plane 

For assigned density profiles and scaled truncation radii, Hj, Ej, two un- 
knowns remain: the fractional mass, rrv (or m), and the fractional radius, y^ 
(or y). The combination of Eqs. fl44a|) and fll03bj) yields: 



w ( ex t) 



8 1 

(V) = g— f {(fi)sd - Cij[(^)mas] 2 } ; (104) 



where the function, w^ ext \ depends on Hj and y\ conformably to Eqs. (I40bl) . 



fTTO . and fl45bl . In the case under discussion, Eq. fll04[) may be conceived 



as a link between the unknowns, rrv and y\ At this stage, one additional 
relation is needed (CV08). 

The mere existence of a fundamental plane (Djorgovski and Davis, 1987; 
Dressier et al., 1987) indicates that structural properties in elliptical galaxies 
span a narrow range, suggesting that some self-regulating mechanism must 
be at work during formation and evolution. In particular, projected light 
profiles from elliptical galaxies exhibit large degree of homogeneity and may 
well be fitted by the r 1//4 de Vaucouleurs law. Accordingly, a narrow range 
may safely be expected also for fractional masses of elliptical galaxies and 
the assumption, m = const, appears to be a viable approximation. This is 
the reason for which the sample used (N = 16) is made of only elliptical 
galaxies, extracted from larger samples (N = 25; iV = 48) of early-type 
galaxies investigated within the SAURON project (SIV; SX). 

Then the fractional radius, y^ (or y), remains as the sole unknown, which 
can be determined by solving Eq. H104[) with numerical techniques. If no 
solution exists, no model galaxy corresponds to the selected sample object 
or, in other words, the related density profiles provide no fit to the data, and 
some input value has to be changed. 

Values of parameters needed for representing model galaxies on the (O^Ty 
^T p ) plane, are listed in Tab.[HJ the galaxy stellar mass, Mf, the galaxy scaling 
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radius, rj; the inclination angle, i, of the best fitting two- integral Jeans model 
(SIV); the anisotropy parameter, 5, determined from the solution of the 
dynamic models, supposed to be axisymmetric (ai = a^) and isotropic on the 
equatorial plane (an = 022) (SX); the intrinsic axis ratio, e, deduced from the 
computed inclination, under the assumption of axisymmetric configurations 
(SX); the dimensionless energy, Cij, defined by Eqs. (I103p . The kinematic 
classification is listed again to get more insight. The galaxy stellar mass 
within the effective radius is calculated as: 

Mi L Mi/L 

L^lO lo m /L ' 

- = exp 10 {-0.4[I T - (to - M) 
and the galaxy effective radius is calculated as: 

Re Re 

kpc arcsec 
Wc = 6XPl ° 

where the factor, 206.265, is related to the choice of measure units. For 
further details refer to an earlier attempt (CV08). 

The values of the reduced variables, fiy, ^T p , ^Tt, are determined via 
Eqs. (I47p and (I50e[) . Both HH and HN macrogases have been considered, for 
the following values of parameters. Scaled truncation radii (both finite or 
infinite): (H^H,) = (ki,kj); ki = 5, 10, 20, +00; kj = 5, 10, 20, +00. Frac- 
tional masses: m = 10, 20. Under the working hypothesis of an analogy 
between VDW gases and macrogases, the (O^Ty^p) plane may be divided 
into three parts: (i) a reversed bell-shaped region where two phases, gas and 
stars, coexist and the lower point coincides with the critical point (hereafter 
quoted as the GS region); (ii) a region limited by the left boundary of the 
reversed bell-shaped region and the rising side of the critical macroisother- 
mal curve, both branching off from the critical point, where only stars are 
present (hereafter quoted as the S region); (iii) a region limited by the right 
boundary of the reversed bell-shaped region and the rising side of the critical 
macroisothermal curve, both branching off from the critical point, and the 
coordinate axes, where only gas is present (hereafter quoted as the G region). 
If the density profiles are only slightly affected in time, the evolution of a 
galaxy on the (O^Tv^p) plane is represented by a track, starting from the G 
region and ending within the GS or the S region. The critical point is the 
sole which is common to the three regions. 



(105a) 

-4.11]} ; (105b) 



d 



206.265 Mpc 
to - M 



(106a) 
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Table 8: Parameters calculated from the data listed in Tab. [7] for a sam- 
ple (N = 16) of elliptical galaxies, extracted from larger samples of early- 
type galaxies investigated within the SAURON project (S IV, N = 25; S X, 
N = 48), which are used in the current paper. Column captions: (1) 
NGC number; (2) galaxy stellar mass deduced from luminosities and mass- 
luminosity ratios (in / band), conformably to Eqs. (1105)) ; (3) galaxy scaling 
radius, calculated using Eqs. ( 11011) and (11061) ; (4) inclination angle of the best 
fitting two- integral Jeans model (S IV); (5) anisotropy parameter, determined 
from the solution of the dynamic models, supposed to be axisymmetric (S X); 
(6) intrinsic axis ratio, deduced from the computed inclination, under the 
assumption of axisymmetric configurations (SX); (7) dimensionless energy, 
defined by Eqs. ( 11031) ; (8) kinematic classification (F - fast rotator; S - slow 
rotator) (SX). For further details refer to the parent papers (S IV; SX) and 
an earlier attempt (Binney, 2005). 
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The position of a model galaxy on the (O^fy^p) plane is affected by er- 
rors of different kind, due to: (1) scatter around mean values listed in Table 
[TJ (2) scatter in fitting observed to model density profiles; (3) uncertainty on 
the determination of the critical point. The third contribution may safely be 
neglected with respect to the other ones. Values related to the first contribu- 
tion are not completely found in literature (to the knowledge of the author). 
The second contribution could be determined using fitting procedures, pro- 
vided light distributions are available and light traces stellar mass in elliptical 
galaxies. In summary, error calculation on the position of sample objects on 
the (O^Tv^p) plane would be cumbersome, and perhaps of little meaning. 

A notable simplification can be attained if model galaxies instead of 
sample objects are considered, in dealing with only errors of the first kind 
mentioned above. For fixed fractional mass, m, the uncertainty on Jj( p = 
(m/m c ) 2 , is negligible with respect to ^Ty = Vc/y, which is determined by 
solving Eq. fll03a|l for assigned density profiles. The combination of Eqs. (llOip . 
(11021) . and (1103bl) yields an expression of the dimensionless energy, Cy, in 
terms of observables listed in Table [7J as: 
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(107a) 
(107b) 
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where some symbols have been changed to gain simplicity. An inspection 
of Eq. ( I107ajl shows that the dimensionless energy, Cy, is monotonically in- 
creasing with increasing £ 4 , ^2, C3; an d decreasing ( 5 , ( 7 ; and vice versa. 
Establishing the trend with the remaining variables, d, i, 5, demands further 
considerations. 

The function, B(e), defined by Eq. (11021) . is monotonically decreasing in 
the domain, < e < 1, where vr = B(0) > B(e) > B{1) = 2. Accord- 
ingly, B(e) is monotonically decreasing with decreasing d =< e± > and/or 
decreasing i, and vice versa. 

With regard to the anisotropy parameter, S, the following identity holds: 



3-5 



1 — 5 cos 2 i 



cos i + 



sin 2 i 



2 cos i 



3-5 



(10* 



and the special inclination angle, io, which makes null the fraction within 
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brackets, is the solution of the equation: 

sin 2 i — 2 cos 2 iq = 3 sin 2 i — 2 = ; 



(109) 



the result is: 



smz = 




i = 0.955 316 6 = 54.735 61° 



(110) 



accordingly, the fraction on the left-hand side of Eq. fllOSp is monotonically 
decreasing for increasing 5 in the range, iq < i < tt/2, and is monotonically 
increasing for increasing S in the range, < i < io, while no dependence on S 
occurs in the special case, i = io- Finally, Eqs. (I103al) . (11081) . and (11091) show 
that the dimensionless energy, c^, is monotonically increasing or decreasing 
for increasing S according if < i < i or i < i < tt/2, respectively, and vice 
versa, while no dependence on 5 occurs in the special case, i = io- 
The above results may be reduced to a single relation, as: 



1 1 C 4 T AC 4 
1163.335 GC5± ACs 

• ex Pl0 [0.4(C 6 T ACe) - 0.2(C 7 ± A( 7 ) - 6.644] 

' arcsin{{l - [1 - (Ci ± ACi)] 2 } 1 ^/ sin(i ± Ai)} 

f 2(C 2T AC 2 ) 2 {3-[^±sgn(^-z )A^]}(C3TAC 3 ) 2 ). 
' \ sin 2 (z ± Ai) l-[6± sgn(z - i )A5] cos 2 (i ± Ai) J ' 



where upper and lower signs correspond to the lower and upper c^- value, 
respectively, and sgn is the sign function, defined as sgn(x) = x/\x\, x ^ 0; 
sgn(0) = 0. It is worth emphasizing that Eq. (II 111) makes an exact for- 
mulation of the uncertainty on the dimensionless energy, c^-, as defined by 
Eq. (I103al) . On the contrary, standard linear and quadratic error propagation 
formulae apply to any kind of functions allowing Taylor series development, 
but are approximate instead of being exact. 

To the knowledge of the author, part of the errors on the right-hand 
side of Eq. (II lip are not available in literature. For the inclination angle, 
i, and the anisotropy parameter, 5, the reason is in that they depend on a 
reference dynamic model (S IV; S X) and cannot be specified. The following 
values are found: AC4/C4 = 0.17 (SIV), 0.20 (SX); ACe/Ce = 0.13 (SIV); 
ACs/Cs — 0.10, strongly dependent on the assumptions made (SIV); A£ 7 = 
0.09 - 0.33, with a value listed for each sample object (SIV). In this view, 
it seems better starting from assigned values of the dimensionless energy 
relative error, Acij/cij, and numerically evaluate the related uncertainty on 
the position of a selected model galaxy on the (O^tv^p) plane, to visualize 
the trend. 
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4.5 Results 



With regard to HH macrogases, model galaxies corresponding to sample ob- 
jects listed in Tables [7] and [8] are represented on the (O^ty^p) plane of Fig.lTBl 
for different choices of scaled truncation radii, Si, Sj, and fractional mass, 
m. The critical macroisothermal curve (left) and the boundary of the GS 
region (right) are also plotted for each case. The critical point is marked 
by a composite symbol. Two sample objects cannot be modelled for low 
inner scaled truncation radii (Sj = 5) and, for this reason, related cases are 
not considered. Lower and upper symbols of the same kind correspond to 
m = 10, 20, respectively. 

Under the working hypothesis of an analogy between VDW gases and 
macrogases, modelled elliptical galaxies are expected to lie in the S region 
or slightly outside the S region within the GS region at most. An inspec- 
tion of Fig. [13] shows the following: (1) model galaxies with low fractional 
mass (m = 10) and/or no truncation radii (S — > +oo) lie below the criti- 
cal macroisothermal curve, in the G region, and for this reason cannot be 
accepted; (2) about one half of model galaxies with low outer scaled radii 
(Sj = 5) lie well inside the GS region, and for this reason cannot be accepted; 
(3) more than one half of model galaxies with larger scaled radii (Sj = 10, 20; 
Sj = 10, 20) lie within the S region, and for this reason are accepted. 

With regard to viable cases, the plot of Fig.[T3]is repeated in Fig.dU where 
model galaxies are distinguished according if their parent sample object is a 
fast (squares) or a slow (diamonds) rotator. The related scaled truncation 
radii (from top to bottom) are (S i; Sj) = (10, 10), (10, 20), (20, 10), (20, 20), 
and the fractional mass is m = 20. The curves are as in Fig.fTBl 

Restricting to viable cases, the plot of Fig. [12] is repeated in Fig.[TS], where 
the effect of assigned errors in dimensionless energy, Acij/cij = 5%, 10%, 
15%, 20%, labelled on each panel, on model galaxies, is represented. More 
specifically, the position of model galaxies is marked by diamonds, and the 
change due to lowered (cy — Acjj) and increased (q, + Acjj) dimensionless 
energy, is marked by Greek and St. Andrew's crosses, respectively. Scaled 
truncation radii are, from top to bottom, (Sj, Sj) = (10, 10), (10,20), (20,10), 
(20,20), and the fractional mass is m = 20 in all cases. In the special cases, 
(Sj, Sj) = (10, 10), (10,20), and Ac^/c^ = 20%, a sample object (NGC 4473) 
cannot be modelled for lowered dimensionless energies, Qj — Acy. In general, 
lowered and increased dimensionless energies, Cjj, make model galaxies shift 
on the left and on the right, respectively, from their position on the (O^Ty^p) 
plane. Accordingly, a fraction of model galaxies enter the G region and the 
GS region, respectively, and the fit could be improved by changing the input 
parameters, Sj, Sj, and m. 
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With regard to HN/NH macrogases, model galaxies corresponding to 
sample objects listed in Tables [7] and [8] are represented on the (O^Ty^p) 
plane of Fig. [16] for different choices of scaled truncation radii, Hi, Hj, and 
fractional mass, m. The critical macroisothermal curve (left) and the bound- 
ary of the GS region (right) are also plotted for each case. The critical 
point is marked by a composite symbol. Two sample objects cannot be 
modelled for low inner scaled truncation radii (Hj = 5) and, for this rea- 
son, related cases are not considered. The same holds, to a larger extent, 
for cases (Hj, Hj) — > (+oo, +oo), due to an infinite mass of the NFW density 
profile. Lower and upper symbols of the same kind correspond torn = 10, 20, 
respectively. 

Under the working hypothesis of an analogy between VDW gases and 
macrogases, modelled elliptical galaxies are expected to lie in the S region 
or slightly outside the S region within the GS region at most. An inspection 
of Fig. [TBI shows the following: (1) model galaxies with low fractional mass 
(m = 10) and/or large outer scaled truncation radius (Hj = 20) lie (at least 
partially) below the critical macroisothermal curve in the G region, and for 
this reason the related cases cannot be accepted; (2) more than one half of 
model galaxies with large fractional mass (m = 20) and/or low outer scaled 
truncation radius (Hj = 5) lie well inside the GS region, and for this reason 
the related cases cannot be accepted; (3) more than one half of model galaxies 
with large fractional mass (m = 20) and outer scaled radius (Hj = 10), or 
low fractional mass (m = 10) and outer scaled radius (Hj = 5) lie within the 
S region, and for this reason the related cases are accepted. 

With regard to viable cases, the plot of Fig.[16]is repeated in Fig.[T71 where 
model galaxies are distinguished according if their parent sample object is a 
fast (squares) or a slow (diamonds) rotator. The related scaled truncation 
radii (from top to bottom) are (Hj,Hj) = (10, 10), (10, 5), (20, 10), and the 
fractional mass is m = 20 for Hj = 10, and m = 10 for Hj = 5. The curves 
are as in Fig.fTBI 

Restricting to viable cases, the plot of Fig. [16] is repeated in Fig. [181 where 
the effect of assigned errors in dimensionless energy, Acij/cij = 5%, 10%, 
15%, 20%, labelled on each panel, on model galaxies, is represented. More 
specifically, the position of model galaxies is marked by diamonds, and the 
change due to lowered (cjj — Acjj) and increased (Qj + Acjj) dimensionless 
energy, is marked by Greek and St. Andrew's crosses, respectively. Scaled 
truncation radii and fractional masses are (Hj,Hj,m) = (10,10,20), (10,5,10), 
(20,10,20), from top to bottom. 

In the special case, (Hj, Hj, m) = (10, 5, 10), and Acij/cij = 20%, a sample 
object (NGC 4486) is out of scale on the right for increased dimensionless 
energies, Cjj + AQj. In the special cases, (Hj, Hj, m) = (10, 5, 10), (10,10,20), 
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and Acij/cij = 20%, a sample object (NGC 4473) cannot be modelled for 
lowered dimensionless energies, — Ac^-. In general, lowered and increased 
dimensionless energies, c^, make model galaxies shift on the left and on the 
right, respectively, from their position on the (O^fv-^p) plane. Accordingly, a 
fraction of model galaxies enter the G region and the GS region, respectively, 
and the fit could be improved by changing the input parameters, Sj, Ej, and 
m. 

4.6 Discussion 

Current cosmological models imply large-scale celestial bodies, such as galax- 
ies and clusters of galaxies, are embedded within nonbaryonic dark haloes, 
where the two subsystems interact only via gravitation. Accordingly, large- 
scale celestial bodies can be modelled as macrogases where macrovolume, 
macropressure, and macrotemperature can be defined, and a counterpart of 
the VDW theory for ordinary gases can be developed. Sufficiently steep den- 
sity profiles, as in HH and HN/NH macrogases, show a similar trend with 
respect to VDW gases: the macroisothermal curves are nonmonotonic with 
two extremum points (one maximum and one minimum) above a critical 
macrotemperature, and are monotonic below a critical macrotemperature. 
The critical macroisothermal curve is characterized by a single extremum 
point, where the maximum and the minimum coincide yielding a horisontal 
inflexion point. On the other hand, sufficiently mild density profiles, such as 
in UU and PP (CV08) macrogases, show only nonmonotonic macroisother- 
mal curves with two extremum points (one maximum and one minimum) 
where no critical macroisothermal curve exists. 

A generic macrogas equation of state is formulated in terms of dimension- 
less variables normalized to critical values (or conveniently chosen in absence 
of the critical point). Similar to what has been done in dealing with the 
reduced VDW equation for ordinary gases (e.g., LL67, Chap. VIII, §85), the 
states of two large-scale celestial bodies with equal JfCy, ^f p , ^Tt, or y, m, 0, 
for assigned scaled truncation radii, Sj, Ej, and belonging to the same family 
of macrogases, may be defined as corresponding states. The mere existence 
of a macrogas equation of state yields the following result. 

Law of corresponding states. Given two large-scale celestial bodies be- 
longing to the same family of macrogases with assigned scaled truncation 
radii, S i; Ej, the equality between two among three reduced variables, fiv, 
^t p , ^Tt, or y, m, <p, implies the equality between the remaining related re- 
duced variables i.e. the two macrogases are in corresponding states. 

The law of corresponding states cannot be extended to macrogases with 
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different scaled truncation radii, as shown in Figs.E] and [T21 A possible 
explanation may be the following. Contrary to ordinary gases, bounded by 
rigid walls which have no influence on the equation of state, macrogases are 
confined by "gravitational" walls appearing in the equation of state via the 
potential energy terms which, in turn, depend on the scaled truncation radii. 

Ordinary gases exhibit monotonic isothermal curves where the central 
part of the related VDW isothermal curve, including the extremum points, 
is replaced by a horisontal line and a phase transition occurs therein. With 
regard to macrogases, the existence of a phase transition and monotonic 
macroisothermal curves of the kind considered, must necessarily be assumed 
as a working hypothesis by analogy with VDW gases. The phase transition 
must be conceived between gas and stars, and the (O^Tv-^p) plane may be 
divided into three parts, namely (i) the G region, where only gas exists; (ii) 
the S region, where only stars exist; (iii) the GS region, where both gas and 
stars exist. In this view, model elliptical galaxies are expected to lie within 
the S region or slightly outside the boundary between the S and the GS 
region at most. 

It is the case for different models related to both HH and HN/NH macro- 
gases, where acceptable values of scaled truncation radii and fractional masses 
are used, as shown in Figs. [13] and [16j The assumption of universal fractional 
mass for sample objects is not a limit of the model, which equally holds 
assigning different fractional masses to different sample objects. 

It can be seen from Figs.[T3] and [T7| that fast rotators lie within the S 
region, while slow rotators are close (from both sides) to the boundary be- 
tween the S and the GS region. This dichotomy could be interpreted by the 
different nature of the two classes of sample objects. More specifically, fast 
rotators seem consistent with elliptical galaxies with disky isophotes, which 
experienced minor mergers and accreted a significant amount of gas (SX), 
suddenly turned into stars. On the other hand, systematically more mas- 
sive (with the exception of NGC 4458 and, marginally, NGC 3608, see Table 
[8]) slow rotators may be related to elliptical galaxies with boxy isophotes, 
which experienced major gas-rich mergers, or sequences of mergers, and reg- 
ulation by the feedback of a powerful central active galactic nucleus (SX), in 
some cases allowing gas survival as e.g., diffuse, still undetected interstellar 
medium. 

Accordingly, in the (O^Ty^p) plane fast rotators and a fraction of slow 
rotators are expected to lie within the S region, and the remaining part 
of slow rotators to be placed in the GS region, as shown in Figs.dH and 
[T71 The following trend is also exhibited: fast rotators are systematically 
on the left with respect to slow rotators, with the exception of NGC 4458, 
which is the sole sample slow rotator with low mass, possibly due to having 
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experienced minor instead of major mergers. In this view, NGC 4458 should 
be considered as a peculiar fast rotator, where low rotation could be due to 
special configurations e.g., still undetected counter- rotating components, as 
observed in the disk of NGC 4450 (e.g., SX). 

The reduced variable, ^Ty = Xy/Xy c: v i a Eqs. ( 142b"|) and (I50ej) is pro- 
portional to the fractional truncation radius, 1/y = Ri/Rj. The above men- 
tioned dichotomy, exhibited by fast and slow rotators in the (O^Ty^p) plane, 
implies a larger fractional truncation radius, y = Rj/Ri, for fast rotators with 
respect to slow rotators. This result could be interpreted as due to different 
formation mechanisms: minor mergers would produce larger contraction of 
the baryonic matter, while major mergers would make (possibly via active 
galactic nuclei) smaller contraction of the baryonic matter, yielding a larger 
or smaller fractional radius, respectively. 

5 Conclusion 

In the current attempt, two-component large-scale celestial bodies where 
the subsystems interact only via gravitation, are conceived as macrogases 
bounded by "gravitational" walls. The macrogas equation of state is formu- 
lated in terms of macrovolume, macropressure, and macrotemperature, which 
are dimensionless variables. For sufficiently steep density profiles, which fit 
to observed elliptical galaxies (or more generally, spheroid components) and 
to simulated nonbaryonic dark matter haloes, macroisothermal curves on 
the [OXyX v ) plane show an analogy with VDW isothermal curves exhibited 
by VDW gases. More specifically, a critical macroisothermal curve exists, 
below and above which the macroisothermal curves are monotonic and non- 
monotonic (with two extremum points, one maximum and one minimum), 
respectively. The critical macroisothermal curve is characterized by a single 
extremum point (a horisontal inflexion point), which is the critical point. 

Contrary to ordinary gases, macrogases cannot be tested in laboratory, 
and for this reason a working hypothesis is inescapable. By analogy with 
ordinary gases, real macroisothermal curves are supposed to occur instead 
of their theoretical counterparts (deduced from the macrogas equation of 
state), where the central part containing the extremum points is replaced by 
a horisontal line along which a phase transition takes place. The intersection 
between a selected theoretical macroisothermal curve and its real counter- 
part, yields two regions of equal area. The phase transition is assumed to be 
gas-stars instead of vapour-liquid as in ordinary gases. Accordingly, the first 
quadrant of the (O^Tv^p) plane is divided into three parts, namely (i) the 
G region, where only gas exists; (ii) the S region, where only stars exist; (iii) 
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the GS region, where both gas and stars exist. 

For selected density profiles and scaled truncation radii, Ej, Ej, the macro- 
gas equation of state depends on three parameters, JjCy, $p, $t, or the 
fractional truncation radius, y, the fractional mass, m, and the fractional 
energy, 0. If elliptical galaxies and their hosting nonbaryonic dark haloes 
are conceived as macrogases, a selected model is accepted only if a whole 
set of sample objects (from which radii, masses, and rsm velocities can be 
determined) lies within the S region or slightly outside the boundary be- 
tween the S and the GS region at most. The sample used (CV08, N = 16) 
is extracted from larger samples of early-type galaxies investigated within 
the SAURON project (SIV, iV = 25; SX, iV = 48). The position of model 
galaxies on the (O^Ty^f p ) plane is determined through the following steps: 
(i) select SAURON data of interest; (ii) calculate the parameters appearing 
in the virial equations; (iii) make a correspondance between model galaxies 
and sample objects; (iv) represent model galaxies on the (O^ty^fp) plane. 

The main results found in the present investigation may be summarized 
as follows. 

(1) A new numerical algorithm has been used for determining the critical 

point of selected HH and HN/NH macrogases, improving earlier re- 
sults (CV08). In particular, the critical point exists for all the cases 
considered. 

(2) A principle of corresponding states rigorously holds for selected den- 

sity profiles and scaled truncation radii, and to a first extent only for 
selected density profiles. 

(3) The following models (on a total of 20) can be accepted in the above 

mentioned sense: (S^S^m) = (10,10,20), (10,20,20), (20,10,20), 
(20, 20, 20), with regard to HH macrogases, and (Sj, Sj, m) = (10, 5, 10), 
(10, 10,20), (20, 10,20), for HN/NH macrogases. The values of model 
parameters may be changed by the occurrence of systematic errors. 

(4) Fast rotators exhibit larger fractional truncation radii with respect to 

slow rotators, which makes the former lie within the S region and the 
latter close (from both sides) to the boundary between the S and the 
GS region, with regard to acceptable models. This dichotomy could 
be interpreted in terms of a different evolution related to fast and slow 
rotators, where gas is currently absent in the former and can be present 
(even if still undetected) in the latter. 
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Appendix 

A Tidal potential energy profile factors 

The tidal potential energy for homeoidally striated ellipsoids related to simi- 
lar and similarly placed boundaries, depends on the reference fractional mass, 
m\ the fractional scaling radius, y\ and the functions, w^ nt \r]), z// ext ) (77) , 
r] = Ei/y^ = ^j/y, expressed by Eqs. ( 14"5"|) . Let the two subsystems be de- 
noted as A and B regardless of what is the outer and what is the inner. 
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Conformingly, Eqs. (|42|) read: 

E B DBA E A E B t r B R-B . , 

— = -r- ; ^ba = -r- = — ; Vba = ~t ' 2/ba = ^- > i ; 

"A yj^ ^ Vba r A Ka 

t t M B M B , s 

U = VbaZb ; m B A = — T ; mBA = XT ; ( 112a ) 

£a 2/ab Sb | i?^ 

— = n - ; vab = | = — ; = -r ; vab = > l ; 

as yJi B 2Mb rjj its 

ft, t M A M A ni9M 

£b = 2/ab^a ; m AB = — r ; ttiab = -j-r- ; (112b) 



B 



according if i? B > i?^ or i?^ > respectively. 

In dealing with sequences of configurations where the scaled truncation 
radii, E A and E B , are kept unchanged, the combination of Eqs. (I112al) and 
ffmbl) yields: 

VbaVab = VbaVab ; (H3) 
for any pair of configurations belonging to opposite sides of the sequence, 
with respect to x/ba = Vab = 1- 

In the special case of equal scaled density profiles, Fa = F B = F, and 
equal scaled truncation radii, Ea = E B = E, Eqs. fj45l) reduce to: 

^(E,yU = ['^ F(0^^^ ; y BA = Vba > 1 (H4a) 

(ext) (2,yL) = jf^ FCyLO-^pe^ ; Vba = Vba>1 (H4b) 

(int H-,^B) = f Mfl ^(0^^^e ; y\ B = y A B>l ;(H4c) 

(ext) (s,</L) = jf^^y^)-^^ ; y ] AB = yAB>i (H4d) 



re 



w 



in the special case where the fractional scaling radius coincides for both 
configurations, y BA = Vabi the combination of Eqs. f)114al) and (jll4cP ; (11 14b|) 
and flll4dl> : yields: 

w^\~,y BA )=w^\~^ AB ) ; Vba = Vab \ (H5a) 
w^\E,y BA ) =w^XE,y AB ) ; y BA = y AB ; (115b) 

regardless of the scaled truncation radius, E. 

If, on the other hand, E A ^ E B , the upper integration limits are E A /y BA 
and E B jy\ B for Eqs. ffTT4aj) and flll4bl> : ffTTicj) and flll4dl) : respectively. 
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Then y BA = y AB implies equal integrands but different upper integration lim- 
its, while 2>a/Uba = ^b/Vab implies equal upper integration limits but dif- 
ferent integrands, with regard to Eqs. fl!14ajl and flll4c|) : (1114b|) and (1114dl) ; 



respectively. Accordingly, Eqs. (I115ap and (1115bj) no longer hold in the case 



under discussion, unless Ea — > +00, Eb — > +00, which erases the dependence 
on E A or E B , regardless of the value of \im^ AtSB )^ +00 (E B /E A ) . 

More specifically, Eqs. (1114]) where E = Ea, E b ; yt = y BA , Vab'i reduce 

to: 



■=->+oo Jo af 



s->+oo Jo df 



d£ 
y ] = y>l ; (116a) 

~ v " ^ d£ 
y ] = y>l ; (116b) 



replacing y^ with 1/y^ yields: 

„M = pp^W^H ; = ( 117 a) 

^4)^(3^^ <ii7b) 

which, choosing £/y^ as integration variable, and keeping in mind that inte- 
grals are independent of integration variables, is equivalent to: 

w dnt) fl\ ..t / + °° ,,.t^ dF(0 



^ \¥j = y] I °° F ^^r^ ; y] = y * 1 ; (us*) 

(^)=l/ t J ( +OO F(0^|^ede ; yt = y >i ; (iisb) 



and the combination of Eqs. fll 16aj) and (1118bj) : (111 6b|) and (jllBaP ; produces: 



e nt) f^) =« xt) (2/ t ) ; (H9a) 



wt t) (^]=y j wt t \y^ ; (119b) 



where, on the other hand, 1/y' — 1/y < 1 is outside the domain, and the 
role of the two components should be interchanged therein, according to 
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Eqs. fll!4cj) and (I114dj) . which makes the above result only mathematically 
relevant. 

Turning to the special case, F A = Fb = F, Ea = H# = H, the last 
impl ying {u A ) S e\ = (^B)sei = ^sei, and taking, in addition, = y\ B = 1, 
Eqs. f)114p reduce to: 

w^ t \~,l)=w^ t \~,l)=w(Z,l) ; (120) 

and the fractional virial potential energy, expressed by Eq. (j47aj) . reduces to: 

, mt^-(9/8M5,l) t 
z/ sc i - (9/8)mTty(e., 1) 

^ei = ^ /V(0d£ ; (121b) 



16 Jo 

u>(s,i)= Tf(0^^ ; ( 121c ) 

jo d-4 

F(S) = ; (121d) 

where Eqs. (1121bl) and (1121d[) follow from the definition of ^ se i and F(£), 
respectively. For further details refer to earlier attempts (e.g., Roberts, 1962; 
Caimmi and Secco, 1992; CV08). 

Integrating by parts Eq. (1121 cl) and combining with (1121bl) and (1121dl) . 
yields: 

tw(S,l) = ~iw ; (122) 

finally, substituting Eq. (I122p into (II 2 lap produces: 

= m f ; yt = 1 ; (123) 

which is a general result for subsystems with equal scaled density profiles, 
F A (£) = F B (Z) = F(£),0<£<S. 



B Dimensional macrogas equation of state 

Let macrogases be defined as large-scale collisionless fluids with the following 
properties: (i) particles are identical mass points; (ii) particle number is 
extremely large; (iii) particle motions obey Newton laws of mechanics; (iv) 
particle collisions are absent; (v) particle interactions obey Newton law of 
gravitation. 

Under the assumption of homeoidally striated density profiles, by use of 
Eq. (j4ip . the virial theorem reads: 

2E kin - u scl G ^ 2 B = ; (124) 
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let the macrovolume, Vm, the macropressure, pyi, the macrotemperature, T M 
and the mass weighted rms velocity, o"m, be defined as: 

An 

~3 

GM 2 

Pm 



M 



-CLlCL2ds 



M 



0"M 



" 3iy ^kin , 



(125) 
(126) 
(127) 



V M 



:i28i 



where N is the total number of particles, k the Boltzmann's constant, and 
the index, M, means macrogas. In particular, Eq. (I127p discloses that the 
macrotemperature, Tm, coincides with the temperature of an ideal gas with 
particle number, N, and translational kinetic energy, E^ m . 

Typical values for galaxies are M — 3 • lO n m , N = 3 ■ 10 11 , <tm = 
lOOv^kms -1 , which yields via Eqs. (II27J) and fTl28|) : &T M = Ma M /(3iV) = 
10 14 -1.99-10 33 erg = 1.99 -10 47 erg, and T M = 1.99- 10 47 /(1.38- 10- 16 )K = 1.44- 
10 63 K = 1440 Km, where where Km = 10 60 K is the macrodegree, assumed as 
macrotemperature unit. 

The combination of Eqs. (I59>(|4TD and ffT24l)- ffT28|) yields: 



PmVi 



M 



47T 



GM 
ai<T M 



47T- 



(129) 



NkT M ai<j M 2z/ scl 5 

which may be conceived as a compressibility factor. 

In presence of a similar, similarly placed, homeoidally striated density 
profile, by use of Eqs. (141 p - ( 143]) . the virial theorem for the subsystem under 
consideration reads: 



2(-E«)km — 



{al)i 



(^u)sel 



1 + 



(^•uujtid 
( z/ u)sel 



B 



and the macrogas equation of state is: 



(Pu)m(V u ) 



M — — 



4tt[(i/ u ) 



u 



^«( z/ «)sel-B 

v = j,i ; 



1 + 



l^mjtid 



K) 



sol 



N u k(T u ] 



M 



(130) 

(131a) 
(131b) 



where i and j denote the inner and the outer subsystem, respectively. The 
compressibility factor is: 



(pJm(K)m 
N u k(T u ) M 

u = i,j ; 



47T- 



GM U 



au)i[(^u)u] 



An 
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^«( I/ u)sel-B 
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( z/ «i))tid 
(^u)sel 



132a) 



(132b) 
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in all cases, the effect of the tidal potential is expressed by the sum within 
square brakets. If the two subsystems were infinitely distant one from the 
other, (i/ u „)tid = 0, and the above results reduce to their one-component 
counterparts. 

In this view, one-component macrogases should be conceived as "ideal" 
and two-component macrogases as "VDW", the related equations of state 
resembling ideal and VDW gas equation of state, respectively. In fact, the 
mass ratio, m = Mj/Mi, and the axis ratio, y = (oj-)i/(aj)i = (a j) 2/(^)2 = 
(0^)3/(04)3, appear in the explicit expression of (i/^tid- Owing to Eqs. (1125P 
and (j!26p . the following relations hold: 



{Vj) M _ ) >{",):', _ , 



(V^m (ai)i(ai) 2 (ai) 3 

(Pj)m = (Mj) 2 (Oi)j(Oi) 2 (Oi)3 = 

(pi) M ~ (M^ 2 (a J ) 2 (a i ) 2 (a J ) 3 ~ y A 

N j( T j)M _ (Ej) km _ ^ 

Ni{Ti) M (Ei) 

kin 



y 6 ; (133) 

(134) 
(135) 



which show the profile factor, (^ w )tid> depends on the fractional macrovolume 
and the fractional macropressure, via m and y. 

The combination of Eqs. ffT27j) . ffTMjl . ffT33l) . and ffT34D yields: 
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(136) 



which, using Eqs. fH5l) . may be cast under the equivalent form: 
1 9 E { (uj) maa 1 ^M(X v ,5i) 



^1 KVi\ 



X P X V n -; 3r 7 sV ^ Sel " p " v = E,)X T ; (137a) 



1 - I ph^l 1 X^X^W^iXy,^) 
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v J mas 



(137b) 



(Pj)m _ m 2 _ „ _ (^)m „, v _ iy,-(Tj 



that is the macrogas fractional equation of state. A simpler and more intu- 
itive choice, adopted in the text, is X p = m 2 , Xy — 1/y, but the connection 
with the fractional macropressure and fractional macrovolume is lost in this 
case. 
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Figure 2: Same as in Fig.[TJ (right panel), where the occurrence (within the 
bell-shaped area bounded by the dashed curve) of saturated vapour is consid- 
ered. Above the critical isothermal curve (T = T c ) the trend is similar with 
respect to ideal gases. Below the critical isothermal curve and on the right 
of the dashed curve, the gas still behaves as an ideal gas. Below the critical 
isothermal curve and on the left of the dashed curve, the liquid shows little 
change in volume as the pressure rises. Within the bell-shaped area bounded 
by the dashed curve, the liquid phase is in equilibrium with the saturated 
vapour phase. A reduced volume implies smaller saturated vapour fraction 
and larger liquid fraction at constant pressure, and vice versa. The VDW 
equation of state is no longer valid in this region. The dashed curve (includ- 
ing the central branch) is the locus of intersection between VDW and real 
isothermal curves, the latter being related to constant pressure where liquid 
and vapour phases coexist. The dotted curve is the locus of VDW isothermal 
extremum points. 
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Figure 3: A specific (T/T c = 0.85) VDW and corresponding real isothermal 
curve. The above mentioned curves coincide within the range, V < Va 
and V > Ve- The VDW isothermal curve exhibits two extremum points: a 
minimum, B, and a maximum, D, while the real isothermal curve is flat within 
the range, Va < V < Ve. Configurations related to the VDW isothermal 
curve within the range, Va < V < Vb (due to tension forces acting on 
the particles yielding superheated liquid), and Vb < V < Ve (due to the 
occurrence of undercooled vapour), may be obtained under special conditions, 
while configurations within the range, Vb < V < Vb, are always unstable. 
The volumes, Va and Ve, correspond to the maximum value in presence of 
the sole liquid phase and the minimum value in presence of the sole vapour 
phase, respectively. The regions, ABC and CDE, have equal area. For further 
details refer to the text. 
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Figure 4: Macroisothermal curves related to IUU (left panel) and AUU 
(right panel) macrogases, respectively. Macroisothermal curves (from bot- 
tom to top) correspond to X T = 0.85, 0.90, 0.95, 1.00, 1.05, 1.10. No crit- 
ical macroisothermal curve exists, above or below or above which the ex- 
tremum points disappear. The coordinates, Xy = X v , X p = X^, X^ = Xp, 
may be conceived as normalized to their fictitious critical counterparts, 

X\r = X\r = 1, X n = X\ = 1, Xt = Xrp = 1. 

Vc Vc ' Pc p c I J c ± c 
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Figure 5: Same as in Fig.H] (right panel). The loci of the extremum points are 
represented as dotted lines. The loci of the intersections between actual and 
real macroisothermal curves are represented as dashed lines. The absence of 
the critical macroisothermal curve makes a band-like instead of a bell-shaped 
region exist on the plane. The coordinates, Xy = Xy, X p = XL Xt = X^, 
may be conceived as normalized to their fictitious critical counterparts, Xy c = 
Xy c = 1) X Pc = X+ c = 1, X Tc = X^ c = 1. 
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Figure 6: A specific (X T = 0.85) AUU and corresponding RUU macroisother- 
mal curve. The above mentioned curves coincide within the range, X v < 
X Va , Xy > Xy E . The AUU macroisothermal curve exhibits two extremum 
points: a minimum, B, and a maximum, D, while the RUU macroisothermal 
curve is flat within the range, Xy A < Xy < Xy E . The regions, ABC and 
CDE, have equal area. The coordinates, X v = X Y , X p = Xj, X T = X^, 
may be conceived as normalized to their fictitious critical counterparts, 
x v c = x Vc = 1, X Pc = X^ c = 1, X Tc = X^ c = I. 
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Figure 7: Macroisothermal curves (X p = X p /X Pc vs. Xy — Xy/X Vc ) related 
to IHH (left panels) and AHH (right panels) macrogases, respectively, for 
different values of scaled truncation radii, (S^S^), labelled on each panel. 
Macroisothermal curves (from bottom to top) correspond to Xt = Xt/Xt c = 
0.90, 0.95, 1.00, 1.05, 1.10, 1.15. The limit, (S^E,-) -> (+oo,+oo), makes 
only little changes. Owing to Eqs. f l42aj) and fl5T|) . X v = Xvy^/y and X^ = 
^(m^/m) 2 , which implies Xy = Xy and X^ = X p . 
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Figure 8: AHH macroisothermal curves (^T p = X p /X Pc vs. ^Ty — Xy/X Vc ) 
for different choices of scaled truncation radii, (Sj, Sj), labelled on each panel. 
Macroisothermal curves (from bottom to top) correspond to ^Tt = X^/Xt c = 
0.90, 0.95, 1.00, 1.05, 1.10, 1.15. The limit, (S^E,-) -> (+oo,+oo), makes 
only little changes. RHH macroisothermal curves, when different from their 
AHH counterparts, lie within the larger reversed bell-shaped region in each 
panel. The loci of intersections between AHH and RHH macroisothermal 
curves are represented as trifid curves, where the left branch corresponds 
to ^V A , the right branch to ^fy E , and the middle branch to $v c - The crit- 
ical point is the common origin. The loci of AHH macroisothermal curve 
extremum points are represented as dotted curves starting from the critical 
point, where the left branch corresponds to minimum points and the right 
branch to maximum points. Owing to Eqs. (I42a|) and (IBTj) . X v = Xyy^/y 
and X^ = X v {m) j 'm) 2 , which implies ^T v = JfCy and ^Tp" = ^f p . 
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Figure 9: AHH macroisothermal curves for scaled truncation radii, (Hj,Hj), 
as in Fig. [HI with regard to the special case, ^Tt = X^/X Tc = 1.15. The limit, 
(Hi, Hj) —> (+00, +00), makes only little changes, as shown by comparison 
with the dashed curve. Owing to Eqs. (jl2aj) and X\ r = Xyyt/y and 

XjJ = X p (m"'"/m) 2 , which implies JjCy = JfCy and ^ = ^t p . 
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Figure 10: Macroisothermal curves (^T p = X p /X Pc vs. ^fy = Xy/XyJ) 
related to IHN/INH (left panels) and AHN/ANH (right panels) macrogases, 
respectively, for different values of scaled truncation radii, (H i; Ej), labelled 
on each panel. Macroisothermal curves (from bottom to top) correspond to 
X T = X T /X Tc = 0.90, 0.95, 1.00, 1.05, 1.10, 1.15. The limit, (E h Ej) 
(+oo, +oo), makes only little changes. Owing to Eqs. (142 al) and f[5"Tj) . X v = 
Xvy^/y and X^ = X^m^ /m) 2 , which implies Jf\ = Xy and JK^ = ^T p . 
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Figure 11: AHN/ANH macroisothermal curves (^T p = X p /X Pc vs. ^Cy = 
Xy/Xv c ) for different choices of scaled truncation radii, Ej), labelled 
on each panel. Macroisothermal curves (from bottom to top) correspond to 
X T = X T /X Tc = 0.90, 0.95, 1.00, 1.05, 1.10, 1.15. The limit, (3*5,-) 
(+oo,+oo), makes only little changes. RHN/RNH macroisothermal curves, 
when different from their AHN/ANH counterparts, lie within the larger re- 
versed bell-shaped region in each panel. The loci of intersections between 
AHN/ANH and RHN/RNH macroisothermal curves are represented as trifid 
curves, where the left branch corresponds to ^V A , the right branch to ^ty E , 
and the middle branch to j[y c . The critical point is the common origin. The 
loci of AHN/ANH macroisothermal curve extremum points are represented 
as dotted curves starting from the critical point, where the left branch corre- 
sponds to minimum points and the right branch to maximum points. Owing 
to Eqs. (I42ap and fl5T|) . X y = Xyy^/y and X^ = X p (m t /m) 2 , which implies 

-^v = -^v and ^ft = JfC p . 
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Figure 12: AHN/ANH macroisothermal curves (^T p = X p /X Pc vs. ^Ty = 
Xy/Xy c ) for scaled truncation radii, (Hj, Ej), as in Fig. [HI with regard to 
the special case, JCt = Xt/Xt c = 1.15. The limit, (Sj,Sj) — > (+oo,+oo), 
makes only little changes, as shown by comparison with the dashed curve. 
Owing to Eqs. (H2aj) and (15T1). X\ T = X v y^/y and = X p (mVm) 2 , which 
implies ^f v = ^Ty and t = ^T p . 
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Figure 13: Elliptical galaxies listed in Tables [7] and modelled as HH macro- 
gases for different choices of scaled truncation radii, Sj, Sj, and fractional 
mass, m. The critical macroisothermal curve (left) and the boundary of the 
GS region (right) are also plotted for each case. Symbol caption and line 
style: (H^H,) = (10,5) - crosses, full; (10,10) - asterisks, dotted; (10,20) 
- diamonds, dashed; (20,5) - triangles, dot-dashed; (20,10) - squares, long- 
short-dashed; (20,20) - St. Andrew's crosses, long-dashed; (+oo, +oo) - dots, 
full. Lower and upper symbols of the same kind are related to m = 10, 20, 
respectively. The composite symbol marks the critical point. Cases where 
Hj = 5 make two galaxies unable to be modelled, and for this reason are not 
considered. 
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Figure 14: Same as in Fig.[T3]for scaled truncation radii, (H*, Hj) = (10, 10), 
(10, 20), (20, 10), (20, 20), from top to bottom, and fractional mass, m = 
20, where model galaxies are distinguished according if their parent sample 
object is classified as fast (squares) or slow (diamonds) rotator. 



68 



0.0 



0.2 



0.-1 



0.6 



X»/X,„ 



0.0 



0.2 



0.-1 



0.6 



X,/X„ 




Figure 15: Same as in Fig. [12] for different choices of scaled truncation 
radii, (H*, Ej) = (10,10), (10,20), (20,10), (20,20), from top to bottom, 
and fractional masses, m = 20, where dimensionless energies, cy, defined 
by Eq. f)103al) . are lowered to cy — Acy (Greek crosses) and increased to 
Cij + Acij (St. Andrew's crosses) with respect to their original values (dia- 
monds), by a factor equal to 5%, 10%, 15%, 20%, respectively. In the last 
case, a sample object (NGC 4473) cannot be modelled for lowered values, 
Clj - A Cij , and {E^Ej) = (10, 10), (10,20). 
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Figure 16: Elliptical galaxies listed in Tables [7] and [HI modelled as HN/NH 
macrogases for different choices of scaled truncation radii, Sj, Sj, and frac- 
tional mass, m. The critical macroisothermal curve (left) and the boundary 
of the GS region (right) are also plotted for each case. Symbol caption 
and line style: (Sj, Ej) = (10,5) - crosses, full; (10,10) - asterisks, dotted; 
(10,20) - diamonds, dashed; (20,5) - triangles, dot-dashed; (20,10) - squares, 
long-short-dashed; (20,20) - St. Andrew's crosses, long-dashed. Lower and 
upper symbols of the same kind are related to m = 10, 20, respectively. The 
composite symbol marks the critical point. Cases where Ej = 5 make two 
galaxies unable to be modelled, and for this reason are not considered. The 
same holds, to a larger extent, for cases (Sj,Sj) — > (+oo, +oo), due to an 
infinite mass of the NFW density profile. 
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Figure 17: Same as in Fig.[16]for scaled truncation radii, (H*, Hj) = (10, 10), 
(10, 5), (20, 10), from top to bottom, and fractional mass, m = 20, Ej = 10, 
and m = 10, Sj = 5, where model galaxies are distinguished according if 
their parent sample object is classified as fast (squares) or slow (diamonds) 
rotator. 
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Figure 18: Same as in Fig. [TH] for different choices of scaled truncation radii 
and fractional masses, (Hj,Hj,m) = (10,10,20), (10,5,10), (20,10,20), from 
top to bottom, where dimensionless energies, c^-, defined by Eq. f)103al) . are 
lowered to Cjj — Acy (Greek crosses) and increased to Cy+Acjj (St. Andrew's 
crosses) with respect to their original values (diamonds), by a factor equal 
to 5%, 10%, 15%, 20%, respectively. In the last sample object (NGC 

4473) cannot be modelled for lowered values, q,j — Acy, and (H^H^m) = 
(10, 5, 10), (10,10,20), and a sample object (NGC 4486) is out of scale on the 
right for increased values, c^- + Acy, and (Hi, Ej, m) = (10, 5, 10). 
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